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EDGE STATES IN HONEYCOMB STRUCTURES 
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Abstract. An edge state is a time-harmonic solution of a conservative wave 
system, e.g. Schrodinger, Maxwell, which is propagating (plane-wave-like) 
parallel to, and localized transverse to, a line-defect or “edge”. Topologically 
protected edge states are edge states which are stable against spatially local¬ 
ized (even strong) deformations of the edge. First studied in the context of the 
quantum Hall effect, protected edge states have attracted huge interest due to 
their role in the field of topological insulators. Theoretical understanding of 
topological protection has mainly come from discrete (tight-binding) models 
and direct numerical simulation. In this paper we consider a rich family of 
continuum PDE models for which we rigorously study regimes where topolog¬ 
ically protected edge states exist. 

Our model is a class of Schrodinger operators on with a background 
two-dimensional honeycomb potential perturbed by an “edge-potential”. The 
edge potential is a domain-wall interpolation, transverse to a prescribed “ra¬ 
tional” edge, between two distinct periodic structures. General conditions are 
given for the bifurcation of a branch of topologically protected edge states 
from Dirac points of the background honeycomb structure. The bifurcation is 
seeded by the zero mode of a one-dimensional effective Dirac operator. A key 
condition is a spectral no-fold condition for the prescribed edge. We then use 
this result to prove the existence of topologically protected edge states along 
zigzag edges of certain honeycomb structures. Our results are consistent with 
the physics literature and appear to be the first rigorous results on the exis¬ 
tence of topologically protected edge states for continuum 2D PDE systems 
describing waves in a non-trivial periodic medium. We also show that the 
family of Hamiltonians we study contains cases where zigzag edge states exist, 
but which are not topologically protected. 


1. Introduction and Outline 

This paper is motivated by a remarkable physical observation. When two dis¬ 
tinct 2-dimensional materials with favorable crystalline structures are joined along 
an edge, there exist propagating modes, e.g. electronic or photonic, whose energy 
remains localized in a neighborhood of the edge without spreading into the “bulk”. 
Furthermore, these modes and their properties persist in the presence of arbitrary 
local, even large, perturbations of the edge. An understanding of such “protected 
edge states” in periodic structures has so far mainly been obtained by analyzing 
discrete “tight-binding” models and from numerical simulations. In this paper we 
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prove that edge states arise from the Schrodinger equation for a class of potentials 
that have many features (not all) in common with the relevant experiments. A cen¬ 
tral role is played by a spectral “no-fold” condition. In the case of small amplitude 
(low-contrast) honeycomb potentials, this reduces to a sign condition of a particu¬ 
lar Fourier coefficient of the potential. A combination of numerical simulation and 
heuristic argument suggests that if the “no-fold” condition fails, then edge states 
need not be topologically protected. Let us explain these ideas in more detail. 

Wave transport in periodic structures with honeycomb symmetry has been an 
area of intense activity catalyzed by the study of graphene, a single atomic layer 
two-dimensional honeycomb structure of carbon atoms. The remarkable electronic 
properties exhibited by graphene 16 23 32 46 have inspired the study of waves in 
general honeycomb structures or “artificial graphene” in electronic 37 and photonic 
contexts. One such property, observed in electronic and photonic 
systems with honeycomb symmetry is the existence of topologically protected edge 
states. Edge states are modes which are (i) pseudo-periodic (plane-wave-like or 
propagating) parallel to a line-defect, and (ii) localized transverse to the line-defect; 
see Figure Topological protection refers to the persistence of these modes and 
their properties, even when the line-defect is subjected to strong local or random 
perturbations. In applications, edge states are of great interest due to their potential 
as robust vehicles for channeling energy. 

The extensive physics literature on topologically robust edge states goes back 
to investigations of the quantum Hall effect; see, for example, 
the rigorous mathematical articles 0[|[^|4^. In 


20 21,41 44 and 


19 34 a proposal for realiz¬ 


ing photonic edge states in periodic electromagnetic structures which exhibit the 
magneto-optic effect was made. In this case, the edge is realized via a domain wall 
across which the Faraday axis is reversed. Since the magneto-optic effect breaks 
time-reversal symmetry, as does the magnetic field in the Hall effect, the resulting 
edge states are unidirectional. 

Other realizations of edges in photonic and electromagnetic systems, e.g. be¬ 
tween periodic dielectric and conducting structures, between periodic structures 
and free-space, have been explored through experiment and numerical simulation; 
see, for example 24 29 35 43 In the context of tight-binding models, the 
existence and robustness of edge states has been related to topological invariants 
(Chern index or Berry / Zak phase E) associated with the “bulk” (infinite periodic 
honeycomb) band-structure. 

We are interested in exploring these phenomena in general energy-conserving 
wave equations in continuous media. We consider the case of the Schrodinger equa¬ 
tion on idtij} = Hip, and study the existence and robustness of edge states 
of time-harmonic form: ip = Our model consists of a honeycomb back¬ 

ground potential, the “bulk” structure, and a perturbing “edge-potential”. The 
edge-potential interpolates between two distinct asymptotic periodic structures, 
via a domain wall which varies transverse to a specified line-defect (“edge”) in 
the direction of some element of the period lattice, Kh- In the context of honey¬ 
comb structures, the most frequently studied edges are the “zigzag” and “armchair” 
edges; see Figure 


Our model of an edge is motivated by the domain-wall construction of 19 34 


A difference is that we break spatial-inversion symmetry, while preserving time- 
reversal symmetry. Hence, the edge states - though topologically robust - may 
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travel in either direction along the edge. In [^[I^ we proved that a one-dimensional 
variant of such edge-potentials gives rise to topologically protected edge states in 
periodic structures with symmetry-induced linear band crossings, the analogue in 
one space dimension of Dirac points (see below). We explore a photonic realization 
of such states in coupled waveguide arrays in [^ . 

Our goal is to clarify the underlying mechanisms for the existence of topologically 
protected edge states. In Theorem |7.3| we give general conditions for a topologically 
protected bifurcation of edge states from Dirac points of the background (bulk) 
honeycomb structure. The bifurcation is seeded by the robust zero mode of a 
one-dimensional effective Dirac equation. A key hypothesis is a spectral no-fold 
condition for the prescribed edge, assumed to be a rational edge. In one-dimensional 
continuum models 12 , this condition is a consequence of monotonicity properties 
of dispersion curves. For continuous d-dimensional structures, with d > 2, the 
spectral no-fold condition may or may not hold; see Sectionj^ Moreover, by varying 
a parameter, such as the lattice scale of a periodic structure, one can continuously 
tune between cases where the condition holds or does not hold; see Appendix [A| 
In Theorem |8.2| and Theorem |8.5| we verify the spectral no-fold condition for the 
zigzag edge, for a family of Hamiltonians with weak (low-contrast) potentials, and 
obtain the existence of zigzag edge states in this setting. 

we study the strong binding regime (deep poten- 


In a forthcoming article II 


tials) for a large class of honeycomb Schrodinger operators. We prove that the two 
lowest energy dispersion surfaces, after a rescaling by the potential well’s depth, 
converge uniformly to those of the celebrated Wallace (1947) tight-binding 
model of graphite. A corollary of this result is that the spectral no-fold condition, 
as stated in the present article, is satisfied for sufficiently deep potentials (high 
contrast) for a very large classes of edge directions in Ah (including the zigzag 
edge). In fact, we believe that the analysis of the present article can be extended 
and together with II will yield the existence of edge states which are localized, 
transverse to arbitrary edge directions t)i S A^. This is work in progress. For a 


detailed discussion of examples and motivating numerical simulations, see 10 


The types of edge states which exist for edges generated by domain walls stand 
in contrast to those which exist in the case of “hard edges”, i. e. edges defined by 
the tight-binding bulk Hamiltonian on one side of an edge with Dirichlet (zero) 
boundary condition imposed on the edge; see parenthetical remark in Figure 
In this case, it is well-known that zigzag (hard) edges support edge states, while 
armchair (hard) edges do not support edge states; see, for example. 


17 


Finally, we believe that failure of the spectral no-fold condition implies that 
there are no topologically protected edge states, although there is evidence that 
there are meta-stable edge states, which are localized near the edge for a long time; 
see Section [TH 


1.1. Detailed discussion of main results. Let Ah = Zvi 0 Zv 2 denote the 
regular (equilateral) triangular lattice and A)( = Zki 0 Zk 2 denote the associated 
dual lattice, with relations k; • = 2Trdim, l,rn = 1, 2. The expressions for k; and 

Vm are displayed in Section [2731 The honeycomb structure, H, is the union of two 
interpenetrating triangular lattices: A 0 A^ and B 0 Ah] see Figures]^ and 

A honeycomb lattice potential, V (x), is a real-valued, smooth function, which 
is Ah— periodic and, relative to some origin of coordinates, inversion symmetric 
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Figure 1. Edge state - propagating (plane-wave like) parallel to a zigzag 
edge (Rvi) and localized transverse to the edge. 


Zigzag edge 




Figure 2. Bulk honeycomb structure, H = (A -h A^) U (B -h A^). Top 
panel: Zigzag edge (blue line), Rvi = {x : k 2 • x = 0}. Shaded region is the 
fundamental domain of the cylinder, corresponding to the zigzag edge. 

Bottom panel: Armchair edge (blue line), R (vi -|- V2) = {x : (ki — k2) • x = 
0}. Fundamental domain of the cylinder, corresponding to the armchair 

edge, also indicated. (Darkened vertices are sites at which zero-boundary 
conditions are imposed in tight-binding models of “hard” edges.) 

















EDGE STATES IN HONEYCOMB STRUCTURES 


5 




Figure 3 . Left panel: A = (0, 0), B = (-^,0). The honeycomb structure, 
H is the union of two interpenetrating sublattices: Aa = A + A^ (blue) 
and Ab = B + Ah (red). The lattice vectors {vi,V 2 } generate Ah- Colors 
designate sublattices; in graphene the atoms occupying Aa— and Ab— sites 
are identical. Right panel: Brillouin zone, Bh-, and dual basis {ki,k2}. K 
and K' are labeled. Other vertices of Bh obtained via application of R, a 
rotation by I'k jZ. 


(even) and invariant under a 27r/3 rotation; see Definition 2.4 A choice of period 
cell is ri/j, the parallelogram in spanned by {vi,V 2 }. 

We begin with the Hamiltonian for the unperturbed honeycomb structure: 

ij(o) = -A +V(x). 

The band structure of the A^— periodic Schrddinger operator, is obtained by 

considering the family of eigenvalue problems, parametrized by k G S/t, the Bril¬ 
louin zone: — E)'^ = 0, tl/(x-|-v) = X G R^, v G At- Equivalently, 

'0(x) = e“®*^ ’''I'(x), satisfies the periodic eigenvalue problem: (k) — E(k)) ijj = 

0 and '(/'(x -I- v) = '(/'(x) for all x S R^ and v G A/j, where = —(V -I- *k)^ -|- 

V(x). For each k G Bh, the spectrum is real and consists of discrete eigenvalues 
Eb(k), b > 1, where Ej(k) < Ej+i(k). The maps k i-)- Eb(k) G R are called the 
dispersion surfaces of The collection of these surfaces constitutes the band 

structure of As k varies over Bh, each map k —>■ Eb{k) is Lipschitz continuous 

and sweeps out a closed interval in R. The union of these intervals is the L^( 
spectrum of E[^^\ A more detailed discussion is presented in Section 

A central role is played by the Dirac points of These are quasi-momentum 

/ energy pairs, (K* , A*), in the band structure of at which neighboring disper¬ 
sion surfaces touch conically at a point 13 23 32 . The existence of Dirac points. 


located at the six vertices of the Brillouin zone, Bh (regular hexagonal dual period 
cell) for generic honeycomb structures was proved in 12 13 ; see also The 


quasi-momenta of Dirac points partition into two equivalence classes; the K— points 
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Figure 4. Lowest three dispersion surfaces k = G Bh -^(k) of 

the band structure of = — A-|-F(x), where V is the honeycomb potential: 
V (x) = 10 (cos{ki • x) -|- cos(k 2 • x) -h cos((ki + k 2 ) • x)). Dirac points occur 
at the intersection of the lower two dispersion surfaces, at the six vertices of 
the Brillouin zone, Bh- 


consisting of K, KK and where i? is a rotation by 27r/3 and K'— points con¬ 
sisting of K' = —K, HK' and The time evolution of a wavepacket, with data 

spectrally localized near a Dirac point, is governed by a massless two-dimensional 
Dirac system [l4] . 

Figure displays the first three dispersion surfaces of for a honeycomb 
potential. The lowest two of these surfaces touch conically at the six vertices of Bh 
(inset). Associated with the Dirac point (K*, A*) is a two-dimensional eigenspace 
of K*— pseudo-periodic states, spanj^i, < 1 ) 2 }: 


nWci>,(x) = n.$,(x), 


X S 


j = 1,2 , where <i>j(x-fv) = e*^*"^$j(x), v G A, 


see Definition 3.1 It is also shown in 13 that a A^—periodic perturbation of y(x). 


which breaks inversion or time-reversal symmetry lifts the eigenvalue degeneracy; a 
(local) gap is opened about the Dirac points and the perturbed dispersion surfaces 
are locally smooth. The perturbation of by an edge potential (see O)) takes 
advantage of this instability of Dirac points with symmetry breaking perturbations. 

To construct our Hamiltonian, perturbed by an edge-potential, we first choose a 
vector 0i G A/i, the period lattice, and consider the line ROi, the “edge”. Choose 
t )2 such that Ah = ZOi © Zt) 2 . Also introduce dual basis vectors, .^1 and ^ 2 , 
satisfying ■ Vm = m = 1, 2; see Section|^for a detailed discussion. The 

choice Oi = Vi (or equivalently V 2 ) is a zigzag edge and the choice Pi = Vi + V 2 is 
an armchair edge; see Figure]^ 

Introduce the perturbed Hamiltonian: 


(1.1) =—A+ V{'k) + 5k{S^2 ■ + d/c(d .^2 • x)lF(x). 
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Here, S is real and will be taken to be sufficiently small, and bb(x) is Ah— periodic 
and odd. The function k, defines a domain wall. We choose k to be sufficiently 
smooth and to satisfy k( 0) = 0 and k{() —>■ ±Koo 7 ^ 0 as C —> ± 00 . Without loss of 
generality, we assume Koo > 0, e.g. k{() = tanh(^). We refer to the line KOi as a 
Pi— edge. 

Note that is invariant under translations parallel to the Pi— edge, x 1—7 
X + Pi, and hence there is a well-defined parallel quasi-momentum, denoted fcii. 


Furthermore, transitions adiabatically between the asymptotic Hamiltonian 


jWix.) as .^2 • X —>■ —00 to the asymptotic Hamiltonian = 
(5KoobF(x) as .^2 ■ X —7 00 . In the case where k changes sign once across 
C = 0, the domain wall modulation of bF(x) realizes a phase-defect across the edge 
(line-defect) KPi. A variant of this construction was used in 12 to insert a phase 
defect between asymptotic dimer periodic potentials. 


Suppose has a Dirac point at (K*,i?*). It is important to note that while 
hw is inversion symmetric, is not. For d 7 ^ 0, does not have Dirac 
points; its dispersion surfaces are locally smooth and for quasi-momenta k such that 
if |k — K*| is sufficiently small, there is an open neighborhood of A* not contained 
in the L^(IR^/A?,)— spectrum of H^\'k). This “spectral gap” about E = may 
however only be local about K* f^. If there is a real open neighborhood of A*, 
not contained in the spectrum of (k) = — (V-|-ik)^ -|-I/±(5 kooVF for all k e 
then E[^^ is said to have a (global) omni-directional spectral gap about E = Ei,. 
We’ll see, in our discussion of the spectral no-fold condition, that it is a “directional 
spectral gap” that plays a key role in the existence of edge states; see Section 1.3 
and Definition l7.ll 

Under suitable hypotheses, we shall construct Pi— edge states of E[^^'>, which are 
spectrally localized near the Dirac point, (K*,U*). These are non-trivial solutions 
Ik, with energies A « A*, of the fc||— eigenvalue problem: 

= E'S, 


( 1 . 2 ) 

(1.3) 

(1.4) 
for fcii 


4'(x-l-Pi) = e*^ii Ik(x) (propagation parallel to KPi), 

|'l'(x)| —7 0 , as • x| —7 00 (localization transverse to MPi), 

: K* - Pi. To formulate the eigenvalue problem in an appropriate Hilbert 


space, we introduce the cylinder E = K^/Zpi. If /(x) satisfies the pseudo-periodic 


boundary condition (1.3|, then /(x)e x jg well-defined on the cylinder E. 

Denote by H®(E), s > 0, the Sobolev spaces of functions defined on E. The pseudo¬ 


periodicity and decay conditions ( |1.3[ )-(1.4) are encoded by requiring 41 S Fl^^^(E), 
for some s > 0, where 


Hf=Him ^ /:/(x)e-*^ 




G Fr"(E) 


Thus we formulate the EVP (|1.2|)-(|1.4|) as: 
(1.5) 


= E'^, ^-GH? (E). 


Remark 1.1 (Symmetry relation among K— and K'— points). Note that if 4'(x) = 
e*k x^^x) is a solution of the eigenvalue problem (1.51, then i/jk = 
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where Z(x-f tJi) = Z(x) and Z(x) —0 as |,^2 -x] —>■ oo, is a propagating edge state 
of the time-dependent Schrddinger equation: = ij(‘^)'!/)(x, t) with parallel 

quasi-momentum k\\ = K • Oi. Since the time-dependent Schrddinger equation has 
the invariance 'ip{x,t) i—> 'ip{x,—t), it follows that 

V:K(x,-t) = = ^K'(x,t). 

Thus '0K' (x, t) is a counterpropagating edge state with parallel quasi-momentum, 
fc|| = K' • Pi = —K • Pi. Due to these symmetry considerations and the equivalence 
of K— points: {K, HK, R^K.}, without loss of generality, we henceforth restrict our 
attention to the Dirac point (K,i5*). 


1.2. Summary of main results. 


1.2.1. General conditions for the existence of topologically protected edge states; 
Theorem \ 7.^ and Corollary \ 7.4\ In Theorem |7.3| we formulate hypotheses on 
the honeycomb potential, V, domain wall function, k(C), and asymptotic periodic 
structure, W{x), which imply the existence of topologically protected Pi— edge 


states, constructed as non-trivial eigenpairs 6 ('1'°, i?”) of (1.51 with fcii = K • Pi, 


defined for all |(5| sufficiently small. This branch of non-trivial states bifurcates from 
the trivial solution branch E = 0, E) aX E = i?*, the energy of the Dirac 

point. Ke y am ong the hypotheses is the spectral no-fold condition, discussed below 
in Section |l.3l At leading order in 6, the edge state, d>'^(x), is a slow modulation 


i:)$_(x) in Hfc 


of the degenerate nullspace of — A* 

( 1 . 6 ) ^'‘^(x) « • x)$+(x) -|- a^,-{SXi2 

( 1 . 7 ) E^ = E, + 0{S^), 0 <| 5 |< 1 , 

where $_(- and $_ are the appropriate linear combinations of <i>i and $ 2 ) defined 


in (4.14). The envelope amplitude-vector, a*(C) = (a*_-|-(C), q:*__(C))”, is a zero- 
energy eigenstate, Pa* = 0, of the one-dimensional Dirac operator (see also ( 6.22| )): 

P = -*|Atj||il2k39^ -H 'dtt«(C)CTi, 


where the Pauli matrices aj are displayed in (1.13). Here Aj G C (see (3.9)) depends 
on the unperturbed honeycomb potential, V, and is non-zero for generic V. The 
constant = (‘I’l, is real and is also generically nonzero. P has 

a spatially localized zero-energy eigenstate for any k{C) having asymptotic limits 
of opposite sign at ±oo. Therefore, the zero-energy eigenstate, which seeds the 
bifurcation, persists for localized perturbations of k(C). In this sense, the bifurcating 
branch of edge states is topologically protected against a class of local perturbations 
of the edge. 

Section gives an account of a formal multiple scale expansion, to any order in 
the small parameter, 5, of a solution to the eigenvalue problem (1.5). The expression 


in (1.6) is the leading order term in this expansion. Our methods can be used to 


prove the validity of the multiple scale expansion, at any hnite order. 

Corollary |7.4| ensures, under the conditions of Theorem |7.3[ the existence of 


edge states, 'I'(x; fey) G 77? (S) for all A:|| in a neighborhood of fey = K • Oi, and by 
symmetry (see Remark|1.4l) for all fcy in a neighborhood of fcy = — K • Oi = K' • Di. 
Thus, by taking a continuous superposition of states given by Corollary |7.4[ one 
obtains states that remain localized about (and dispersing along) the zigzag edge 
for all time. 
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Remark 1.2. A key hypothesis in Theorem |7.3| is a spectral no-fold condition at 
(K, A*) for the Oi— edge of the band-structure of —A-\-V. This (essentially) ensures 
the existence of a (S)— spectral gap containing A* for the perturbed 


Hamiltonian, see Dehnition 7.1 and the discussion in Section 1.3 


1.2.2. Theorem \8.^ ' Existence of topologically protected zigzag edge states. We con¬ 
sider the case of zigzag edges corresponding to the choice t)i = Vi, t )2 = V 2 , and 
.^1 = ki, .^2 = ^ 2 . Recall that = Zvi © Zv 2 . The choice Oi = V 2 would lead 
to equivalent results. 

We consider the zigzag edge state eigenvalue problem 


( 1 . 8 ) 




(see also (1.5)), 


with Hamiltonian 


(1.9) = -A + ey(x) + SK { 6k2 ■ x)W(x) 


+5N(5k2 •x)lT(x). 


Here, e and 6 are chosen to satisfy 

(1.10) 0<|(5|<e2<l. 


There are two cases, which are delineated by the sign of the distinguished Fourier 
coefficient, eVi,i, of the unperturbed (bulk) honeycomb potential, eH(x). Here, 

^ 1.1 ^ V{y) dy, 

is assumed to be non-zero. We designate these cases: 


J- ,-L 


In Appendix|^we give two explicit families of potentials, a superposition of “bump- 
functions” concentrated, respectively, on a triangular lattice, A^“\ and a honeycomb 
structure, H, that can be tuned between these two cases by variation of a lattice 
scale parameter. 

Under the condition eUip > 0 (Case (1)) and (|1.10 ), we verify the spectral 
no-fold condition for the zigzag edge in Theorem |8.2 The existence of zigzag 
edge states (Theorem |8.5[) then follows from Theorem 7.3 and Corollary 7.4 In 


particular, for all e and <5 satisfying (1.10) and for each ki\ near K • vi = 2tt/3, the 


zigzag edge state eigenvalue problem (1.5) has topologically protected edge states 


with energies sweeping out a neighborhood of El, where (KjA®) is a Dirac point. 

Remark 1.3 (Directional versus omnidirectional spectral gaps). While the regime of 
weak potentials, implied by (1.10), would at first seem to be a simplifying assump¬ 


tion, we wish to remark on a subtlety for = —A + eV ± 5kooW (e, S small), 

which arises precisely in this regime. It is well-known that for sufficiently weak 
periodic potentials on d > 2, that there are no spectral gaps; this is related to 
the “Bethe-Sommerfeld conjecture” [4,38 39 . Nevertheless, if eVip > 0, and £ and 
6 are related as in (1.10), then a directional spectral gap, i.e. an (S)— spectral 
gap exists; see Theorem |8.3|and Section [T31 
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Figurej^and Figure|^are illustrative of Cases (1) and (2). The simulations were 
done for the Hamiltonian with e = ±10 and 0 < J < 10: 


+ (5K((5k2 • x)lF(x), K,{C,) = tanh(^), 

( 1 . 11 ) ^ ^ 

y(x) = ^cos(ii-’ki • x), lF(x) = ^(— 1)'^^^ sin(i?Ai • x). 

i=o 1=0 

Here, R is the rotation matrix displayed in (2.6). Figure displays, for 

fixed e, the ^fe||= 27 r/ 3 (^)~ spectra (plotted horizontally) of corresponding 

to a range of (5 values (strength / scale of domain wall -perturbation) for Cases 
(1) eVip > 0 (top panel) and (2) eVi^ < 0 (middle and bottom panels). Figure 
^displays, for these cases, the (S)— spectra (plotted vertically) for a range of 
parallel-quasi-momentum, ku. 


Remark 1.4 (Symmetries of /cy i—>■ E{k^\)). Figurej^exhibits some elementary sym¬ 
metries. Since the boundary condition for the EVP ( |1.8[ ), 'k(x± Vi) = e*^ii'k(x) is 
27r— periodicity in /cy, the mapping fcy i—>■ i?(/cy) is 27r— periodic. Furthermore, in¬ 
variance under complex conjugation, implies symmetry of fcy i—>■ E{k\\) about fcy = 0 
and fcy = TT. 


1.2.3. Non-topologically protected bifurcations of edge states. In Case (2), where 
£Vi,i < 0, Theorem 8.4 implies that the spectral no-fold condition fails and we 
do not obtain a bifurcation from the Dirac point. However, through a combination 
of formal asymptotic analysis and numerical computations, we do find bifurcating 
branches of edge states. These branches do not emanate from Dirac points (the no¬ 
fold condition fails), but rather from a spectral band edge. Moreover, as we discuss 
below, these states are not topologically protected; they may be destroyed by an 
appropriate localized perturbation of the edge. Case (2) (eVi,i < 0) is illustrated 
by Figures [^(middle and bottom panels) and Figure]^ (bottom panel). 

In particular, Dirac points occur at the intersecti on o f the second and third 
spectral bands of = —A ± eV(x) (see Theorem |3.5| ), and the failure of the 

spectral no-fold condition implies that an — spectral gap does not open about 
E = El ioi 5 and small. However, for eVi^ < 0 there is a spectral gap between 
the first and second spectral bands of For the choice of edge-potential 

displayed in (1.11) with e = —10, a family of nontrivial edge states bifurcates, 
for 0 < |(5| sufficiently small, from the upper edge of the first (lowest) 
spectral band into the spectral gap (dotted blue curve); see middle panel of Figure 

A bifurcation of a similar nature is discussed in [^ . 

A formal multiple scale analysis clarifies this latter bifurcation. For k G let 
(E®(k), $^(x; k)) denote the eigenpair associated with a lowest spectral band. In 
[To] , We calculate that the edge state bifurcation is seeded by a discrete eigenvalue 
effective Schrodinger operator: 


( 1 . 12 ) His = - 


1 




2mls dC 


7 ^+Qeff(C;«). where 


1 


''eff 


ij = 1.2 




and (5eff(C; >^'{0 + ^ ~ ^^(0) is a spatially localized effective po¬ 

tential, depending on k(C), and constants a and 6, with 6 > 0, which depend on 
V, W and For the above choice of the zigzag edge-potential (middle panel of 
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Figure 5. spectra, where K ■ vi = Itt, of the Hamiltonian 

(i fTnl l) for the zigzag edge Top panoh Case (1) 0. 

Topologically protected bifurcation of edge states, described by Theor em |8.2| 
(dotted red curve), is seeded by zero-energy mode of a Dirac operator | |6.22[ |. 
The branch of edge states emanates from intersection of first and second bands 
{Bi and B 2 ) at E = for 5 = 0; see discussion in Section |l.2. ^ Middle 
panel: Case (2) eVi^\ < 0 with domain wall function k. Spectral no-fold 
condition does not hold. Bifurcation of zigzag edge states from upper endpoint, 
E = E^, of the first spectral band. This bifurcation is seeded by a bound state 
of a Schrodinger operator l |1.12[ l with effective mass mgff < 0 and effective 
potential Qeff(C) (displayed in the inset) and is not topologically protected; 
see discussion in Section |1.2.3| Bottom panel: Case (2) < 0 with 

domain wall function k,^. Bifurcation from upper endpoint of Bi is destroyed. 
Bound states bifurcate from the lower edges of the first two spectral bands. 


Figure]^, we have < 0 and the effective potential Qlg, displayed in the figure 
inset, induces a bifurcation into the gap above the first band. 

Now, we can construct domain wall functions, for which the corresponding 

has no point eigenvalues in a neighborhood of the right (upper) edge of the 
first spectral band; see bottom panel of Figure If k{() is chosen as above, then 
Qeff(C; (1 ~ )j 0 < 0 < 1, provides a smooth homotopy from a Schrodinger 

Hamiltonian for which there is a bifurcation of edge states with domain wall 
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> 0 




Figure 6. Top panel: L| (S)— spectrum of protected states of 

for the case > 0. Bottom panel: (S)— spectrum of non-protected 

states of for the case £:Fi,i <0. V,W and k, are chosen as in {LU}. 

For each fixed fc||, edge states shown in the top panel > 0) arise due to 

a protected bifurcation from a Dirac point displayed in the top panel of Figure 
Those edge states indicated in the bottom panel (eUi,! < 0) arise via an 
edge bifurcation of the type shown in the middle and bottom panels of Figure 
The band edge energies from which this latter bifurcation takes place is 
well-separated from the energy of the Dirac point which, when eVi^i < 0, lies 
within the overlap of the second and third spectral bands. 


k) to one for which the branch of edge states does not exist with domain wall 

K[,). Therefore, this type of bifurcation is not topologically protected; see 


10 


for a 

more detailed discussion. This contrast between topologically protected states and 
non-protected states is explained and explored numerically, in a one-dimensional 
setting in 


27 


1.3. Remarks on the spectral no-fold condition. The spectral no-fold hy¬ 
pothesis of Theorem |7.3| requires that the dispersion curves obtained by slicing 
the band structure (situated in x with a plane through the Dirac point 
(K,iil*) containing the direction K 2 (dual direction to the Oi— edge) do not fold- 
over and fill out energies arbitrarily near if*. This essentially implies that via 
a small perturbation which breaks inversion symmetry (as we do with = 

—A -I- V(x) -I- Sk{S ^2 ■ x)lV(x) for (5 ^ 0) we open a (S)— spectral gap about 
if*. Figureis illustrative. 

In the first row of plots in Figure we consider whether the spectral no-fold 
condition holds at the Dirac point (K, if*) for the zigzag edge, in the two cases: 
(1) eVi^i > 0 and (2) eVi^i < 0, as well as for the armchair edge. The energy level 
if = if* is indicated with the dotted line. In the left panel we see that for the 
zigzag edge, the spectral no-fold condition holds if eVi^i > 0. In this case, there is 
a topologically protected branch of edge states. In the center panel we see that the 
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Zigzag eV^ ^>0 Zigzag eV^ ^<0 





Armchair eV^ 




- 0.2 0 0.2 
A 


Figure 7. Zigzag and armchair slices at the Dirac point (K, of the 
band structure of —A-\-eV -\-5k,ooW for 5 = 0 (first row) and 5 > 0 (second 
row). Insets indicate zigzag and armchair quasi-momentum segments (one¬ 
dimensional Brillouin zones) parametrized by A, for 0 < A < 1. See discussion 
of Section [Til and Theorem 


spectral no-fold condition fails if eVi^i < 0. Finally, in the right panel we see that 
it also fails for the armchair slice. 

The second row of plots in Figure]^ illustrates that the spectral no-fold condition 
controls whether a full — spectral gap opens when breaking inversion symmetry. 

In particular, for <5 > 0, is no longer inversion symmetric. For eVi^i > 0, 

a spectral gap opens about the Dirac point, between the first and second spectral 
bands (see Theorem 3.5). For the zigzag edge with eVi^i < 0 there is no spectral 


gap about the Dirac point. (Note, however, that there is a spectral gap between the 
first and second spectral bands; see the discussion above in Section 1.2.3 ) Similarly, 
for the armchair edge (right panel) there is no spectral gap for (5 > 0 . 

1.4. Are there meta-stable edge states? Consider the Hamiltonian = 

—Ax + V{x)+6k {6^2 ■ x) W{x) (as in ([13), corresponding to an arbitrary rational 
edge, KOi, i.e. Oi = OiVi -|- 61 V 2 , Oi and bi co-prime integers, as introduced in 



finite order in the small parameter 6. 

But is this formal expansion the expansion of a true edge state? We believe the 
answer is no, if the spectral no-fold condition fails. 

Indeed, from Theorem 


4.2 


we have that any t>i— edge state, 'k G 




is a superposition of Floquet-Bloch modes of = —A + V along the quasimo¬ 
mentum segment: K -|- A.^ 2 , |A| < 1/2. The formal expansion of Section how¬ 
ever is spectrally concentrated on Floquet-Bloch components along this segment. 
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which are near the Dirac point, corresponding to |A| ^ 1. If the spectral no-fold 
condition fails, the expansion does not capture the effect of resonant coupling to 
quasi-momenta along this segment “far from K” (corresponding to A bounded away 
from A = 0 in Figure [^. 

Conjecture: Suppose the spectral no-fold condition fails for the Oi— edge MPi. 
Then, has topologically protected long-lived (meta-stable) edge quasi-modes, 
Ip G ^fc||=K Di,ioc(^); generically has no topologically protected edge states. 


1.5. Outline. 

In Section we review spectral theory for two-dimensional periodic Schrddinger 
operators, introduce the triangular lattice, the honeycomb structure and honeycomb 
lattice potentials. 

In Section we define Dirac points and review the results on the existence of 


Dirac points for generic honeycomb potentials from 13,14 

In Section we introduce the notion of an edge or line defect in a bulk (un¬ 
perturbed) honeycomb structure. Honeycomb structures with edges parallel to a 
period lattice direction, have a translation invariance. Thus, an important tool 
is the Fourier decomposition of states which are (localized) in the unbounded 
direction, transverse to the edge, and propagating (plane-wave like) parallel to the 
edge. 

In Section]^ we introduce our class of Hamiltonians, consisting of a bulk hon¬ 
eycomb potential, perturbed by a general line-defect / tJi— edge potential. 

In Section we give a formal multiple scale construction of edge states to any 
finite order in the small parameter 6. 

In Section we formulate general hypotheses which imply the existence of a 
branch of topologically protected Ui— edge states, bifurcating from the Dirac point. 
The proof uses a Lyapunov-Schmidt reduction strategy, applied to a system for the 
Floquet-Bloch amplitudes which is equivalent to the eigenvalue problem. Such a 
strategy was implemented in a ID setting in 12 . First, the edge-state eigenvalue 


problem is formulated in (quasi-) momentum space as an infinite system for the 
Floquet-Bloch mode amplitudes. We view this system as consisting of two coupled 
subsystems; one is for the quasi-momentum / energy components “near” the Dirac 
point, and the second governs the components which are “far” from the 

Dirac point. We next solve for the far-energy components as a functional of the 
near-energy components and thereby obtain a reduction to a closed system for the 
near-energy components. The construction of this map requires that the spectral 
no-fold condition holds. 

In Section [^we consider the Hamiltonian, introduced in Section]^ in the weak- 
potential (low-contrast) regime and prove the existence of topologically protected 
zigzag edge states, under the condition eVi_i > 0. 

In Appendix we give two families of honeycomb potentials, depending on 
the lattice scale parameter, a, where we can tune between Case (1) eVi^i > 0 and 
Case (2) eVi^ < 0 by continuously varying the lattice scale parameter. 

In a number of places, the proofs of certain assertions are very similar to those 
of corresponding assertions in 12 . In such cases, we do not repeat a variation on 


the proof in 12 , but rather refer to the specific proposition or lemma in 12 


1.6. Notation. 
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(1) Vj, j = 1,2 are basis vectors of the triangular lattice in A/j. k^, £ = 1,2 
are dual basis vectors of Aj!^, which satisfy • Vj = 

(2) For m = (toi, m 2 ) G mk = miki + m 2 k 2 . 

(3) tJi = aiVi +a 2 V 2 G A^, oi, 02 co-prime integers. The Oi— edge is KtJi. Oj, j = 

1,2, is an alternate basis for A/j with corresponding dual basis, = 1,2, 

satisfying ■ \3j = 27ri5^j. 

( 4 ) 3 ^ J^(l) I3I = 

(5) B denotes the Brillouin Zone, associated with A^,, shown in the right panel of 
Figure]^ 

( 6 ) (/, 5 ) = !lg- 

(7) X <y if and only if there exists C > 0 such that x < Cy. x y ii and only if 
X <y and y ^ x. 

(8) LP’^(R) is the space of functions F : K —>■ M such that (1 -|- G L^(K), 

endowed with the norm 


lie’ll 


Lp>‘ 


(l + |f)*/"F 


LP( 


u 

3=0 


F 


LP( 


< 00 , 1 < p < 00 . 


(9) For f,g G the Fourier transform and its inverse are given by 

Him = m = 7^ f e-’^ «/(A)dA, T-\g}{X) = g{X) = f 

(,27!) Jjjd 


The Plancherel relation states: f{x)g{x)dx = (27r)‘^ 

(10) (Jj, j = 1, 2, 3, denote the Pauli matrices, where 


(1.13) 





and (73 



1.7. Acknowledgements. We would like to thank I. Aleiner, A. Millis, J. Liu and 
M. Rechtsman for stimulating discussions. 


2. Floquet-Bloch Theory and Honeycomb Lattice Potentials 
We begin with a review of Floquet-Bloch theory; see, for example, [§[^[2§|^. 

2.1. Fourier analysis on L^(]R/A) and L^(E). Let { 01 , 02 } be a linearly inde¬ 
pendent set in and introduce the 

Lattice: A = ZOi © Z 02 = (miOi -I- m 202 : mi,m 2 GZ}; 

(2.1) Fundamental period cell: H = {0iOi + 62^2 ■ 0 < 0^ < 1, j = 1,2}; 

Dual lattice: A* = © Z ^2 = (mil = mi.^i + 1712^2 ■ mi, m 2 G Z}, 

K, ■ Vj = 2Tr6ij, 1 < i, j < 2; 

Brillouin zone: B, a choice of fundamental dual cell; 

Cylinder: E = IR^/Z0i; 

(2.2) Fundamental domain for E: Hs = (riOi + T 2 O 2 : 0 < ti < l,r 2 G K}. 

We denote by L^(H) and L^(Os) the standard spaces on the domains H and 
Hs, respectively. 

Definition 2.1. [The spaces L^(]R^/A) and L^.] 
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(a) L^(IR^/A) denotes the space of functions which are A— periodic: / G 

L^(IR^/A) if and only if /(x -|- tj) = /(x) for all x G 0 G A and 

fGL^n). 

(b) denotes the space of functions which satisfy a pseudo-periodic 

boundary condition: /(x -|- tj) = e**^ *’/(x) for all x G 0 G A and 

G L^(M^/A). For / and g in fg is in L^(M^/A) and we 
define their inner product by 

{f,9)Ll = [ f{^)9(^)dx- 
“ Jn 


Definition 2.2. [The spaces and 

(a) A^(S) = L^(M^/ZOi) denotes the space of functions, which are periodic 
in the direction of Pi: /(x-|-tii) = /(x), for all x G and such that 
/ G where fii; is the fundamental domain for S; see (2.21. 

(S) = denotes the space of functions: 

( 1 ) which are ki\— pseudo-periodic in the direction Ci: 


(b) 


/(x-ht)i) = e*'=ii/(x), forxGM^ and 

( 2 ) such that which is defined on S, is in 

For / and g in Lf. (S), fg is in Lf{'^) and we define their inner product 

by 

(/,5)l2 = [ /(x)g(x)dx. 

'‘II JQs 

The respective Sobolev spaces iJ®(M^/A), H^, and are de¬ 

fined in a natural way. 

Simplified notational convention: We shall do many calculations requiring us 
to explicitly write out inner products like (/, 5 )i 2 ( 2 ) {f,g)i ^2 ( 2 )- We shall 

write these as f(x)g(x) dx rather than as f{x)g{x) dx. 


If / G L^(K^/A), then it can be expanded in a Fourier series: 
(2.3) 


/(x) = ^ U = 




1 




f{y)dy , mil = wiili -f 1712^2 , 


where jnj denotes the area of the fundamental cell, ff. In Section 4.1 we show 
that, if 5 G L'^{E), then it can be expanded in a Fourier series in Pi • x and Fourier 
transform in P 2 • x: 


g{x) = 2 ^ V / " 


2 ti gn(27r^) = 


1 


y „-inAi y 


|tJi A O2I Jy. 


9{y)dy ■ 


2 . 2 . Floquet-Bloch Theory. Let Q(x) denote a real-valued potential which is 
periodic with respect to A. We shall assume throughout this paper that Q G 
(^ 00 ( 1 ^ 2 /A), although we expect that this condition can be relaxed without much 
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extra work. Introduce the Schrodinger Hamiltonian H = —A + Q(x). For each 
k € we study the Floquet-Bloch eigenvalue problem on 

(2.4) Fr$(x;k) = £;(k)$(x;k), x e 

$(x + tj) = e*’^ “<i)(x; k), Vt) S A. 


An solution of (2.4) is called a Floquet-Bloch state. 


Since the k— pseudo-periodic boundary condition in (2.4) is invariant under 
translations in the dual period lattice, A*, it suffices to restrict our attention to 
k S where B, the Brillouin Zone, is a fundamental cell in k— space. 

An equivalent formulation to (2.4) is obtained by setting $(x; k) = e**^ ’‘p(x; k). 
Then, 

(2.5) i7(k)p(x; k) = £’(k)p(x; k), X G K^, p(x -|-o) = p(x; k), 0 G A, 

where H(k) = —(V -I- ik)^ -|- Q{x) is a self-adjoint operator on L^(IR^/A). The 


eigenvalue problem (2.5), has a discrete set of eigenvalues Ei(k) < £’ 2 (k) < • • • < 
A'b(k) < •••, with L^(IR^/A)— eigenfunctions pf,(x;k), b = 1,2,3,— The maps 
k G H I—>■ Ej(k) are, in general, Lipschitz continuous functions; see, for example. 
Appendix A of [^. For each k G H, the set {pj(x;k)}j>i can be taken to be a 
complete orthonormal basis for L^(IR^/A). 

As k varies over B, £'{,(k) sweeps out a closed real interval. The union over 
6 > 1 of these closed intervals is exactly the L^(IR^)— spectrum of — A -|- H(x): 
spec(iJ) = spec (FI(k)) . Furthermore, the set {<i)(,(x; k)};,>i_kGB is complete 
in 


(^b(-;k),/(-))L2(B2)®b(x;k)c?k = ^ / / 6 (k)$t,(x;k)dk 

b>i b>i 


where the sum converges in the norm. 


2.3. The honeycomb period lattice, A^, and its dual, A)(. Consider Ah = 
Zvi 0 Zv 2 , the equilateral triangular lattice generated by the basis vectors: Vi = 
= (^ 1 — 5 )^; see Figure]^ left panel. The dual lattice A)) = Zki © 
Zk 2 is spanned by the dual basis vectors: ki = q{^, ^)^, k 2 = q{^,—^)'^, where 
<7 = ^, with the biorthonormality relations k^ • Vj = 2Tr5ij. Other useful relations 
are: |vi| = |v 2 | = 1, vi • V 2 = |, |ki| = |k 2 | = 9 and ki • k 2 = — \q^- The Brillouin 
zone, Bh, is a regular hexagon in Denote by K and K' its top and bottom 
vertices (see right panel of Figure given by: K = ^ (ki — k 2 ) , K' = —K = 
^ (k 2 — ki). All six vertices of Bh can be generated by application of the matrix 
R, which rotates a vector in clockwise by 27r/3: 



The vertices of Bh fall into two groups, generated by the action of i? on K and 
K': K— type-points: K, i?K = K + k 2 , i?^K = K — ki, and K'— type-points: 
K', RK' = K' - k2, R^K.' = K' + ki. 

Functions which are periodic on with respect to the lattice Ah may be viewed 
as functions on the torus, K^/A/j. As a fundamental period cell, we choose the 
parallelogram spanned by vi and V 2 , denoted H/,. 
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Remark 2.3 (Symmetry Reduction). Let ($(x;k),iJ(k)) denote a Floquet-Bloch 
eigenpair for the eigenvalue problem (2.4) with quasi-momentum k. Since V is real, 


($(x;k) = <i)(x; k), i?(k)) is a Floquet-Bloch eigenpair for the eigenvalue problem 
with quasi-momentum —k. The above relations among the vertices of Bh and the 
A)(- periodicity of: k i-)- £'(k) and k <i)(x;k) imply that the local character of 
the dispersion surfaces in a neighborhood of any vertex of Bh is determined by its 
character about any other vertex of Bh- 

2.4. Honeycomb potentials. 

Definition 2.4. [Honeycomb potentials] Let V be real-valued and V G 
H is a honeycomb potential if there exists Xq S such that V (x) = H (x — Xq) has 
the following properties: 

(VI) V is Ah— periodic, i.e. V(x -|- v) = V{x) for all x G and v G Ah- 
(V2) V is even or inversion-symmetric, ie. V(—x) = V(x). 

(V3) V is TZ- invariant, ie. 7?.[V](x) = V{R*x) = V(x), where, R* is the 
counter-clockwise rotation matrix by 27r/3, ie. R* = R~^, where R is given 


by (2.6). 


N.B. Throughout this paper, we shall omit the tildes on V and choose coordinates 
with Xq = 0. 

Introduce the mapping R ■- 1? ^ 1? which acts on the indices of the Fourier 
coefficients of V: R{mi,m 2 ) = {—m 2 ,mi—m 2 ) and therefore (to 1 , m 2 ) = {m 2 — 
mi, —mi), and R?{mi,m 2 ) = {mi, m 2 )- Any m ^ 0 lies on an R— orbit of length 
exactly three [^. We say that m and n are in the same equivalence class if m 
and n lie on the same 3— cycle. Let S denote a set consisting of exactly one 
representative from each equivalence class. Honeycomb lattice potentials have the 


following Fourier series characterization 13 : 

Proposition 2.5. Let V(x) denote a honeycomb lattice potential. Then, 

V (x) = tq -p ^ Vm cos(mk • x) -|- cos((i?m)k • x) -|- cos((i?^m)k • x) 

mGS 

where mk = miki -|- TO 2 k 2 and the Vm o-re real. 


3. Dirac Points 

In this section we summarize results of [13| on Dirac points. These are conical 
singularities in the dispersion surfaces of Hy = — A-f V(x), where V is a honeycomb 
lattice potential. 

Let K* denote any vertex of Bh, and recall that is the space of K*— pseudo- 
periodic functions. A key property of honeycomb lattice potentials, V, is that Hy 
and TZ, defined in (V3) of Definition 


2.4 


leave a dense subspace of 




invariant. 
Hy commutes with TZ: 


Furthermore, restricted to this dense subspace of 
[TZ, Hy] = 0. Since TZ has eigenvalues 1, t and t , it is natural to split L 
direct sum: 


Ij- into the 




= 






Here, where a = 1, r, t and r = exp(27rj/3), denote the invariant eigenspaces 

of7^: 
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We next give a precise definition of a Dirac point. 


Definition 3.1. Let ^(x) be a smooth, real-valued, even (inversion symmetric) 
and periodic potential on Denote by Bh, the Brillouin zone. Let K. G Bh- The 
energy / quasi-momentum pair (K,iJ*) G Bh x ^ is called a Dirac point if there 
exists 6* > 1 such that: 

(1) if* is a eigenvalue of Hy of multiplicity two. 

(2) Nullspace^iiy — if*jj = span|<i)i(x), $ 2 (x)|, where G Lk r = 

Tibi) and <i) 2 (x) = (C oZ) [<i)i](x) = $i(-x) G { 71^2 = t ^ 2 ), and 

(‘ho) ‘f’ti)L^(n) “ a,b = 1,2. 

(3) There exist Aj 7 ^ 0, Co > 0, and Floquet-Bloch eigenpairs 

ki-G {^b,+i{x;'k),Et^+i{k)) and k i-a (x; k),(k)), 

and Lipschitz functions ej(k), j = 6*,6* -I- 1, where ej(K) = 0, defined for 
|k-K I < Co such that 

Eb^+i{k)-E, = +|Att| |k-K| (1 + e6.+i(k)), 

(3.1) EbA^)-E, = -IAjjI |k-K| (1 + 66,(k)), 

where |ej(k)| < Clk — K|, j = 6*, 6* -I- 1, for some C > 0. 


In 13 , the authors prove the following 

Proposition 3.2. Suppose conditions 1 and 2 of Definition \3.l\ hold and let |c( m)| ing 5 
denote the sequence of Fourier-coefficients o/$i(x) normalized as in Il3|. 

Define the sum 


(3.2) 


mG5 ^ ^ 


Here, S Gif is defined in [13| . If 0, then condition 3 of Definition 3.1 holds 
(see (3.1)/ 


Therefore Dirac points are found by verifying conditions 1 and 2 of Definition |3.1| 
and the additional (non-degeneracy) condition: Aj / 0. 

Furthermore, Theorem 4.1 of and Theorem 3.2 of [13 imply the following 
local behavior of Floquet-Bloch modes near the Dirac pointTTH 


Corollary 3.3. 

(3.3) 

$b,+i(x;k) = 

(3.4) 

$6,(x;k) = 


1 r A# 

(k-K)(i) -hf(k-K)(2) 

4>i(x) -h 4'2 (x) 

72 L |A#| 

|k-K| 

1 r ^ 

(k-K)(i) -hf(k-K)(2) 

$l(x) - 4>2 (x) 

72 L Ajt 

|k-K| 


$ 


( 1 ) 

b.-l-l 


(x;k). 




where (•; k) = C)(|k - K|) in H^{Flh) as |k - K| ^ 0. 


In the next section we discuss the result of 13 , that —A -|- eV has Dirac points 
for generic e. 


^The factor // 


in ||3.3[l-l[3l4|l corrects a typographical error in equation (3.13) of 
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3.1. Dirac points of — A+el^ (x), e generic. The strategy used in to produce 
Dirac points is based on a bifurcation theory for the operator — A + ey(x) acting in 
L^, from the e = 0 limit. We describe the setup here, since we shall make detailed 
use of it. 

Consider —A acting on We note that = |Kp is an eigenvalue with mul¬ 
tiplicity three, since the three vertices of the regular hexagon, Bh- K, RK and i?^K 
are equidistant from the origin. The corresponding three-dimensional eigenspace 
has an orthonormal basis consisting of the functions: ‘&cr(x) = e*^'^pcr(x) G 
g., cr = 1,T, r, defined by 


(3.5) 

We note that 


(cr = l,r,r) 


1 


1 


+ cre*^ 


iK-: 


»k2X , 


1 + aE'^^'^ + ae 


(3.6) 


— {PcnP^) L^(m.‘2/Ah) 


= s. 


In Theorem 5.1 of [13| , the authors proved that for real, small and non-zero e 
and under the assumption that V satisfies the non-degeneracy condition: 


(3.7) 


Di,i = 




[ D(y) dy ^ 0, 

h I J 


that the multiplicity three eigenvalue, = |Kp, splits into 

(A) a multiplicity two eigenvalue, A®, with two-dimensional r®^K t~ eigenspace 
structure, and 

(B) a simple eigenvalue, A®, with one-dimensional eigenspace, a subspace of i- 
For all e sufficiently small, the quasi-momentum pairs (K,ifJ) are Dirac points 

in the sense of Definition 3.1 Furthermore, a continuation argument is then used to 


extend this result from the regime of sufficiently small e to the regime of arbitrary 
e outside of a possible discrete set; see 


13 


and the refinement concerning the 
possible exceptional set of e values in Appendix D of . We first state the result 
for arbitrarily large and generic e, and then the more refined picture for |e| > 0 and 
sufficiently small. 


Theorem 3.4. [Generic e] Let D(x) be a honeycomb lattice potential and consider 
the parameter family of Schrodinger operators: 

= -A-heD(x), 

where V satisfies the non-degeneracy condition ( |3.7| ). Then, there exists Sq > 0, 
such that for all real and nonzero e, outside of a possible discrete subset of M.\ 
{—eo,£o), has Dirac points (K,iff) in the sense of Definition 

Specifically, for all such e, there exists 5* > 1 such that A* = E[ (K) = 
is a K— pseudo-periodic eigenvalue of multiplicity two where 

(1) (a) El is an eigenvalue of of multiplicity one, with correspond¬ 

ing eigenfunction, $^(x). 


3.1 
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(h) El is an eigenvalue of of multiplicity one, with correspond¬ 
ing eigenfunction, $|(x) = x). 

(c) El is not an eigenvalue of H^^'>. 

(2) There exist 5e >0, Cg > 0 and Floquet-Bloch eigenpairs: ($®(x; k), i!^|(k)) 
and Lipschitz continuous functions, Cjfk), j = 6*,^* + 1, defined for |k — 
K| < Jg, such that 

El+,ik)-E%K) = +|Af| |k-K| (1 + e?^+i(k)) and 

(3.8) Elk)-E^K) = - |Af| |k-K| (l + elk)), 
and where 

(3.9) H c{m^El,ef ( • ) • (k + mk) ^ 0 

is given in terms of {c{iLn,El,e)}in^s, the Fourier coefficients of 

<i)f(x;K). Furthermore, |e|(k)| < Cs\k — K|, j = &*,6* + 1. Thus, in a 
neighborhood of the point {k, E) = (K, El) G the dispersion surface is 
closely approximated by a circular cone. 

3.2. Dirac points of —A + ei^(x), e small. In this section we collect explicit 

information on Dirac points for the weak potential regime. 


Theorem 3.5. [Small e] There exists Eq > 0 , such that for all e G = {—eo,eo)\ 
{0} the following holds: 

(1) For e G leo, — A + eT (x) has 

(a) a multiplicity two L^- eigenvalue Ef, where ker{—A -G eV) C Tk ^ 0 

and 

(b) a multiplicity one L^- eigenvalue El, where ker(—A 0 eV) C i- 

(2) The maps e ^ El and e i—>■ El are well defined for all e in the deleted 
neighborhood of zero, 1^^. They are constructed via perturbation theory of a 
simple eigenvalue in ^ and in , respectively. Therefore, El and El 
are real-analytic functions of e G leo ■ Moreover, they have the expansions: 

( 3 . 10 ) El =\K\'^+£{Vop-Vi,i)+0{e^), 

(3.11) El = |K|2 + eCDo.o 0 2Vi,i) 0 0{e‘^). 


(3) If eVi^i > 0, then conical intersections occur between the I'** and 2"'^ dis¬ 


persion surfaces at the vertices of Bh- Specifically, (3.8) holds with 6* = 1. 


(4) If eVi^i < 0, then conical intersections occur between the 2"'^ and dis¬ 
persion surfaces at the vertices of Bh- Specifically, (3.8) holds with 6* = 2. 


For e G I, 


Eq 7 


(3.12) 


= 47r|D?i| 0 0{e) = 47r|vi A V 2 I 0 0{e). 


(6.25) 


The expansions (3.10), (3.11) and (3.12) are displayed in equations (6.22), 
and (6.30) of 13 . 

The intersections of the first two dispersion surfaces for eVi^i > 0, and of the second 
and third dispersion surfaces for eTip < 0, are illustrated in the first two panels of 
Figure along a dispersion slice corresponding to the zigzag edge. 
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4. Edges and dual slices 

Edge states are solutions of an eigenvalue equation on which are spatially 
localized transverse to a line-defect or “edge” and propagating (plane-wave like 
or pseudo-periodic) parallel to the edge. Recall that Kh = Zvi 0 Zv 2 and = 
Zki 0 Zk 2 . We consider edges which are lines of the form ]R(aiVi 0 a 2 V 2 ), where 
(ai,&i) = 1 , i.e. ai and &i are relatively prime. 

We fix an edge by choosing Oi = aiVi 0 I 11 V 2 , where (ai,&i) = 1. Since ai,&i 
are relatively prime, there exists a relatively prime pair of integers: 02,62 such that 
0162 — 02^1 = 1- Set t )2 = a 2 Vi 062 V 2 . It follows that Zt 3 i 0 Zt )2 = Zvi 0 Zv 2 = A/i. 
Since 01&2 — a 2 &i = 1, we have dual lattice vectors ^ 1,^2 € A)), given by 

.^1 = 62ki — a 2 k 2 , ^2 = + 01^27 


which satisfy 

■ t)£' = 2'k5£^£', 1 < £,(' <2. 

Note that Z.fti 0 Z ^2 = ^ki 0 '^^2 = A)(. 

Fix an edge, MOi. In our construction of edge states, an important role is played 
by the “quasi-momentum slice” of the band structure through the Dirac point and 
“dual” to the given edge. 

Definition 4.1. For the edge KOi, the band structure slice at quasi-momentum K, 
dual to the edge Mpi, is defined to be the locus given by the union of curves: 

A ^ £;b(K 0 AJ^ 2 ), |A|< 1 / 2 , &> 1 . 

We give two examples: 

(1) Zigzag: Pi = vi, 02 = V 2 and .^1 = ki and ^2 = k 2 . 

In this case, we shall refer to the zigzag slice. 

(2) Armchair: Di = vi 0 V 2 , ti 2 = V 2 and .^1 = ki and .^2 = k 2 — ki. 

In this case, we shall refer to the armchair slice. 

Figure]^ (top row) displays three cases, for —A 0 eV, where E is a honeycomb 
lattice potential. Shown are the curves A i-A- £'b(K 0 A.^ 2 ), b = I, 2, 3 for (i) eEi.i > 0 
and the zigzag slice (left panel), (ii) eEiq < 0 and the zigzag slice (middle panel), 
and (iii) the armchair slice (right panel). As discussed in the introduction, of 
these three examples, case (i) is the one for which the spectral no-fold condition of 
Definition O holds. 

4.1. Completeness of Floquet-Bloch modes on L^(E). For Oi G Kh, introduce 
the cylinder E = K^/ZOi. Consider the family of states $f,(x;K 0 A.^ 2 )! b > 1 
for A G [0,1] (or equivalently |A| < 1/2) corresponding to quasi-momenta along a 
line segment within Bh connecting K to K 0 .^ 2 - Since .^2 ■ = 0 , all along this 

segment we have K • Di— pseudo-periodicity: 

$f,(x 0 ni;K 0 Ai? 2 ) = e*^^+^^=’ ‘’^$h(x;K 0 Ail 2 ) =e*^ ”i$ 6 (x;K 0 Ail 2 ) ■ 

The main result of this subsection is that any / G is a superposition 

of these modes. 


Theorem 4.2. Let / G k-oi (S) = Then, 
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(1) f can be represented as a superposition of Floquet-Bloch modes of—A + V 
with quasimomenta in B located on the segment = K +A ^2 : ^ 5 j" ■ 


/(x) \ /b(A)‘i>b(x;K + Ai? 2 )dA 

b>i 

1 

(4.1) = e*^-^=’'/6(A)pb(x;K +All 2 )dA, where 

L'-.. 1 ^ i 


b>l “2 


MX) = ($b(-,K + Al^ 2 ),/(-))L? , (E)- 

fe||=K-DiV 


Here, the sum representing e“^^'^/(x)j in (4.1) converges in the 
norm. 

(2) In the special case where V = 0; 


/(x) = E 


(K+mli)-: 


/m(A). 




dX . 


Proof of Theorem 14 We introduce the parameterizations of the fundamental pe¬ 
riod cell n of 14 (x): 

xGll: X = TitJi-I-r 2 t) 2 , 0<ri,r2<l, k* • x = 27rri, 

(4.2) dxi dx2 = |t)iAti2| dri dT2 = |n| dri c?T2 ; 

and of the cylinder S = 

X G E : X = TitJi -I- r 2 t) 2 , 0 < Ti < 1, s G K, • x = 27rri, .^2 ' x = 27rr2, 

(4.3) dxi dx2 = |t)i A 02 ! dri dT2 = |n| dri dT2. 

Let / G L^ii^j^.p^(E) be such that g{-x.) = e~^^'^f{x) is defined and smooth 
on E, and rapidly decreasing. It suffices to prove the result for such /, and then 
pass to all Lfc||=K oi(^) standard arguments. The function g{x) has the Fourier 
representation 

(4.4) = Ott In 


5(x) = 271 V /" 5„(27r^)e*^^" ’"d^e 


271 g„(27r^) = 


1 

tJi A 0 




2 \ JE 


9{y)dy ■ 


The relation (4.4) is obtained by noting that G(ri,r 2 ) = giriPi + T 2 t> 2 ) is 1— 
periodic in ti and in L^(M; dT 2 ), and applying the standard Fourier representations. 
Introduce the Gelfand-Bloch transform 

(4.5) 5(x;A) = 2^ E (27r(m2 + A)) |A| < 1/2 . 

(mi ,m2)GZ2 


Note that x >-)• (ji(x; A) is Ah— periodic and A i-)- g(x; A) is 1— periodic. Using ( |4.5[ ) 
and ( 

(4.6) 


and (4.4), it is straightforward to check that 

1 

5 (x)= F^^--^gMX) dX. 
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Remark 4.3. For any fixed |A| < 1/2, the mapping x i-A g(x; A) is A/j— periodic. We 
wish to expand x i-)- g(x; A) in terms of a basis for Lp‘{Vl), where denotes our choice 
of period cell (parallelogram) for M^/A with A = ZtJi 0 Zt) 2 ; see (2.11. Now the 
eigenvalue problem iJ(k)p^ = on ft with periodic boundary conditions has 

a discrete sequence of eigenvalues, £'j^(k), j > 1 with corresponding eigenfunctions 
P5'(x;k), J > 1, which can be taken to be a complete orthonormal sequence. Recall 
p)^'‘(x;k), b > 1, with corresponding eigenvalues, Fl&(k), the complete set of eigen¬ 
functions of H{\i) with periodic boundary con ditio ns on the elementary period 
parallelogram spanned by {vi, V 2 }; see Section 2.2 By periodicity, (x; k), 6 > 1, 
(initially defined on klh) and p^(x;k), j > 1, (initially defined on 11) can be ex¬ 
tended to all as periodic functions. We continue to denote these extensions by: 
p^(x;k) and p^'*(x;k), respectively. Since A = ZOi © Zt )2 = Zvi © Zv 2 = A/i, 
both sequences of eigenfunctions are Ah— periodic. Thus, in a natural way, we 
can take p)/(x;k) = At p|^'*(x;k), 6 > 1, where Ah is a normalization constant. 
Abusing notation, we henceforth drop the explicit dependence on SI, and simply 
write pt,(x; k) for p^(x; k). 


In view of Remark |4.3| we expand ^(x; A) in terms of the states {p{,(-;K + 

Ai^2)}, b>l: 

(4.7) p(x;A) =^(pb(-;K + Ail2),5(-,A))^2(f2) Pb(x; K + AAl 2 )- 

b>l 

Recall that /(x) = e^^'^g{x). We claim (and prove below) that 

(4.8) (pb(-;K + Al%2),5(-,A))^2(f,) = ($b(-;K + AAl2),/(-))L| • 


The assertions of Theorem 4.2 then follow from (4.61, (4.7) and the claim (4.8): 

1 

/(x) = e*K-5(x) = 




T 

6>1 ’ 


g*Aii2x~(x;A)dA 

/_ 1 
2 

(Pb(-; K + XA2),9{; X))^^) K + X^2)dX 


= E / ' K + ^-*^ 2 ), /(•))l 2 (s) b(x; K + XR,)dX 


b>l * 


where, in the final line we have used that ${,(y; A) = p&(y; Therefore, 

it remains to prove claim (4.8). We shall employ the one-dimensional Poisson 
summation formula: 


(4.9) 


271 ^ / (27r(n + A)) = '^f{y + . 


n^'L 


In the following calculation we use the abbreviated notation: 


P 6 (y;A) =Pb(y;K + AI^2) and $f,(y; A) = $t,(y; K + Al^2)- 
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Substituting (4.4|-(4.5| into the left hand side of (4.8) and applying (4.9) gives 
A),3(-,A)>^2(o) = A)g{y, A)dy 

= |0| / [ + r 202 ; A)27r (27r(m2 + A)) (| |4.2| |) 

= |0| / f Pi,(Tit)i + r2D2; A) ^ 

Jo Jo „ 


TTLl 


Stt ^ (27r(m2 + A)) 


m2 GZ 


e2’r™lDrf.ridT2 
e2’r™iDrf.ridT2 (ra) 


= |0| r r pUtiVi + r2D2; A) ^ ^ (r2 + 

Jo Jo miez m 2 ez 

= |0| f f Pi,(Tit)i + st)2; A)e“^’"®^'’ ^ finds (r2 + m 2 = s) 

■Jo -'K migz 

= |0| / f P 5 (tiBi + st) 2 ; A)e“^'^®^'’G(Ti, s)fiTifis (G(Ti,r 2 ) = g(TiOi +T 2 P 2 ) and l |4.4[ |) 

Jo Jr - 

= J^Pb{y,^)e-^^^^-^giy)dy (by |^). 

Finally, recalling that g = and Ph(y; = $f,(y; A), we obtain 

that 

(P5(-; A),ff(-,A))^ 2 (n) = /(y) dy 

= / ^’h(y;A)/(y) dy = ($b(-;A),/(-))i 2 . 

JE fc||=K.DiV ) 

This completes the proof of claim (4.8) and part 1 for the case where / is smooth 
and rapidly decreasing. Passing to arbitrary / £ is standard. In the case 
where P = 0, the Schrodinger operator reduces to the Laplacian —A. In this case 
the Floquet-Bloch coefhcients of / £ ^fc||=K oi simply its Fourier coefficients: 

/(A) = /(A). Part 2 therefore follows from part 1, completing the proof of Theorem 

10 □ 


Sobolev regularity can be measured in terms of the Floquet-Bloch coefficients. 
Indeed, as in Lemma 2.1 in 12 , by the 2D— Weyf law ^’^(K -|- A.ft 2 ) ~ b for all 
A £ [—1/2,1/2], 6 ^ 1, we have 

Corollary 4.4. Lfc||=K oi(^) ^nd ld/||=K Bi(^)’ s £ N, norms can be expressed 
in terms of the Floquet-Bloch coefficients fbW, b > 1. For f £ Lk^^='K-oi('^) ~ 

iL=K.B,(KVz«i)-- 




fcii =K-Ui (^) 


E 

6>1 * 


l/5(A)pdA , 


E(i+^)' 


|/5(A)pdA . 


b>l 


4.2. Expansion of k i-A i?b(k) along a quasi-momentum slice. Let (K, E*) 

denote a Dirac point as in Definition |3.1| In a neighborhood of the Dirac point, 
the eigenvalues £’f,^(k) and Eb^+i(k) are Lipschitz continuous functions and the 
corresponding normalized eigenmodes, 4)b^(x;k) and 4>b^+i(x;k) are discontinuous 
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functions of k; see 14 . Note however that what is relevant to our construction of 


til— edge states are Floquet-Bloch modes along the quasi-momentum line K-|-A,^ 2 j 
I A| < 1/2. The following proposition gives a smooth parametrization of these modes 
along this quasi-momentum line. 


Proposition 4.5. Let (K, E^) denote a Dirac point in the sense of Definition 
Let {<i>i(x), $ 2 (x)} denote the basis of the L^ = nullspace of 


3.1 


Hv — Ei,I in Definition \3.1\ Introduce the A/j— periodic functions 

(4.10) Pi(x) = e-*^’^4>i(x), P2(x)=e-*^'"<i>2(x). 


For each\X\ < 1/2, there exist eigenpairs X), E±{X)), real analytic 

in X, such that ($a(-; A), ${,(•; A)) = Sab and 


span {$_(x; A), <i>+(x; A)} = span (x; K-|- A.ft 2 ), 4)b,-|-i(x; K-|- A.ft 2 )} • 


Introduce Ah— periodic functions p±(x; A) by 

(4.11) 4>±(x;A) = p±(x; A), {Pa{-; >^)) = Sab, a, 6 G{-f,-}. 


There is a constant (/g > 0 such that for |A| < the following holds: 

(1) The mapping X i—>■ E±{X) is real analytic in X with expansion 

(4.12) E±{X)=E,±\X^\ |J^ 2 | A + P 2 .±(A)A" , 

where Aj G C is given by ( |3.2[ ), |p 2 ,±(A)| < C with C a positive constant 
independent of X. 

(2) Let ^2 = I 32 I = |•^ 2 |• The Ah— periodic functions, p±(x;A), 

can be chosen to depend real analytically on X and so t/ia<|^ 

(4.13) p±(x; A) = P±(x) -h (p±(x,A) G Ll^{]Sfi/Ah), 


where p± (x; 0) = P± (x) is given by 


and (x; 0 ) = (x) is given by 


(4.14) 


$±(x) 


J_ r 32 

1/2 L lAjjl I32I 


4>i(x) ± $ 2 (x) 


Finally, X 1 —>■ <p±(x; A) are real analytic satisfying the bound |clx<p± (x;A)| < 
C'A for all x G Ah, where H = (^ 1 ,^ 2 ), |H| < 2. 


N.B. We wish to point out that the subscripts ± have a different meaning here 
than in [^[^. In [^[^, P±(k) denote ordered eigenvalues, P-(k) < P-(-(k) 
(Lipschitz continuous) with corresponding eigenstates (x; k) (discontinuous at 
k = K) ; see Definition |3.1| and Corollary |3.3[ In Proposition |4.5| and throughout this 
paper P±(A) and 4)±(x; A) refer to smooth parametrizations in A of Floquet-Bloch 
eigenvalues and eigenfunctions of the spectral bands, which intersect at energy P*. 


The factor of in | |4.13[ l corrects a typographical error in equation (3.13) of 


14 
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Proof of Proposition\4.5[ We present a proof along the lines of Theorem 3.2 in 14 ; 


see also |15[|22| . The k— pseudo-periodic Floquet-Bloch modes can be expressed in 
the form <l>(x; k) = k), where p(x; k) is periodic. For k = K -f A.^ 2 , 

consider the family of eigenvalue problems, parametrized by |A| < 1/2: 

(4.15) 7Fv(K -I- Ai^ 2 ) p(x; A) = if(A)p(x;A), 

(4.16) p(x-b v; A) = p(x; A), for all v G , 

where i/y(k) = — (Vx + + l^(x). Degenerate perturbation theory of the 

double eigenvalue if* of iFy(K), yields eigenvalues: if±(A) = if* -b if£^(A), where 

(4.17) ^^±^(A) = ± lAjI |i^ 2 | A+ 0(A"); see [^. 

Denote by Q± the projection onto the orthogonal complement of span{Pi,P 2 }- 
Then, 

iiK(P*) = (iiy(K) Q^L'^{RyAh)^Q±L'^{M.yAh) 

is bounded. Furthermore, via Lyapunov-Schmidt reduction analysis of the periodic 


eigenvalue problem (4.15)-(4.16) we obtain, corresponding to the eigenvalues (4.17), 


the Ah— periodic eigenstates: 

p±(x;A) = (J + i?K(£;*)QT(27Al^2-(V + iK)))x 

(a(A) Pi(x) -b /3(A) P 2 (x)) -b 
Here, the pair a(A),/3(A) satisfies the homogeneous system: 


|2^i 




= 0 , where 


/? 


M{E^^\X) = 


£;(!)+0(A2) _AjA32+0(A2) 

-A#A^ + 0(A2) F;(i)+0(A2) 


see [^. For if^^^ = E^p{X), j = ±, normalized solutions, pj(k;A), j = ±, are 
obtained by choosing: 


o+(A) 

P+W 


( 


A 


J_ _32_ 


75 HF UfT ^ 
0(A) 


1 

■x/5 


o_{A) 

/3_(A) 


/ 


75 l^itl I 32 I 


■0(A) 


V 


1 

’ a/2 


+ C(A) y 


Finally, we note that M.{E^^\X) is analytic in the parameter A. Therefore the 
eigenvalues E^\x) and eigenvectors (q!±(A),/ 3±(A))^ are analytic functions of A; 
see, for example, [^22 . It follows that if±(A) and p±(x;A) are bounded, real 
analytic functions of A G M. This completes the proof of Proposition |4.5[ □ 

5. Model of a honeycomb structure with an edge 


Let V (x) denote a honeycomb potential in the sense of Definition 2.4 In this sec¬ 
tion we introduce a model of an edge in a honeycomb structure. A one-dimensional 
variant of this model was introduced and studied in [9 12,27. 

Let W G (7°°(M^) be real-valued and satisfy the following properties: 


(Wl) W is Ah— periodic, i.e. IT(x-b v) = IF(x) for all x G and v G Ah- 
(W2) W is odd, i.e. IT(—x) = —IF(x). 
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(W3) ) ^ 0, with $1 as in Definition 


3.1 


The non-degeneracy condition (W3) arises in the multiple scale perturbation theory 
of Section m 

Our model of a honeycomb structure with an edge is a smooth and slow interpola¬ 
tion between the Schrodinger Hamiltonians = — Ax -I- H(x) — 6kooW{x.) and 

= —Ax -I- H(x) -I- (5koo1T(x), which is transverse to a lattice direction, say 
t)i. Here, Koo is a positive constant. This interpolation is effected by a domain wall 
function. 


Definition 5.1. We call k(C) € a domain wall function if k(C) tends to 

±Koo as C “> ±oo, and Ti(C) = k^(C) ~ ^ 2 (C) = ^^(C) satisfy: 

(5.1) /(I-|-|C|)“|T^(iC)|d^ < CX3 for some a > 5/2 and f < oo, ^=1,2. 

«/ IR J W. 

Without loss of generality, we assume Koo > 0. 


Remark 5.2. The technical hypotheses in (5.11 are required for the boundedness of 

See also Section 6.6 of 


wave operators used in the proof of Proposition 7.15 


and, in particular, the application of Theorem 6.15. 


12 


Our model of a honeycomb structure with an edge is the domain-wall modulated 
Hamiltonian: 

= — Ax -I- H(x) -I- Sk { 6^2 ■ x) W (x) , 

where /c(C) is a domain wall function. Suppose k(C) has a single zero at C = 0. The 
“edge” is then given by KOi = {x : .^2 • x = 0}. 

We shall seek solutions of the eigenvalue problem 

(5.2) = 

(5.3) 'k(x -I- Ui) = e**^ *’^'I>(x) (propagation parallel to the edge, KUi), 

(5.4) 'k(x) —^ 0 as |x-.ft 2 |— >-00 (localization tranverse to the edge, KtJi). 

In the next section we present a formal asymptotic expansion of Di— edge states 
and in Section we formulate a rigorous theory. 


6. Multiple scales and effective Dirac equations 


We re-express the eigenvalue problem ( |5.2[ )-(5.4 ) in terms of an unknown function 
4' = 4'(x, ^), depending on fast (x) and slow (C = 6^2 ■ x) spatial scales: 


( 6 . 1 ) 


-(Vx + 51^2 ac)VH(x) 


4- 


(5k(C)W(x)4' = 


(6.2) 4'(x-b0i,C) =e*^"i4'(x,C), and 4'(x,C)^0 as (^± 00 . 


We seek a solution to (|6.1|)-(|6.2|) in the form: 

(6.3) 

(6.4) 


E& = i;(o) + ^^(1) j2^(2) ^ ^ 

¥ = (x, C) + (x, C) + 521 /( 2 ) (x, ^) 


The conditions (5.3), (5.4) are encoded by requiring, for i > 0, that 
(x -f t), •) = (x, •) V 0 e A/j, 

C^i/«(x,C)eL2(M^). 
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Substituting the expansions (6.3)-(6.4| in ( |6.1[ ) yields 

[- (Ax + 2S ^ 2 - Vx 5c + 5^2 + ...) + (F(x) + 6 k{0W{x)) 

- + <52^(2) + . . ^^(0) + ^^(1) + J2^(2) ^ ^ q_ 

Equating terms of equal order in j > 0, yields a hierarchy of equations governing 

At order (5° we have that satisfy 


(6.5) 


(^-Ax + E(x) - = 0, 

i/>(°)(x + t), •) = •) Vue A^. 


Equation (6.5) may be solved in terms of the orthonormal basis of the L^(n)— 
nullspace of Ely — E^, in Definition 3.1 namely {$ 1 , $ 2 }. Expansion (4.13) in Propo¬ 


sition 


4.5 suggests that a particularly natural orthonormal basis of the L^[VL)— 


nullspace of Ely — A* is given by {4)+, <!)_}, where 


( 6 , 6 ) 




MB 


Here Aj is given in (3.2), 32 = -^ 2 ^^ and I 32 I = |.^ 2 |- We therefore solve (6.5) 

with 

(6.7) V'^0)(x,C)=a+(C)<i>+(x) + a-(C)4>-(x). 

Proceeding to order we hnd that satisfies 

(-Ax + P(x) - E^) V^(i)(x, C) = g(i)(x, C; 

-I- 0 , •) = •) Vue Ah, 

where 

G(^)(x,C;' 0^°^) = G(^)(x,C;a+,Q!_) 

= 25^0!+ ^2 ■ Vx4>+ -I- 2dt^a- ^2 ■ — «((C)bE(x) (Q!+<i>+ -|- «_$_). 


( 6 . 8 ) 


Viewed as an equation in x, (6.8) is solvable if and only if its right hand side is 
L^(r2;dx)— orthogonal to the nullspace of Hy — A*. This is expressible in terms 
of the two orthogonality conditions: 

= 2 ($j, 1^2 • Vx$+) 5 ca+ + 2 ($j, 1^2 • Vx$-) d^a. 

(6.9) -«(C) [($„H/$+)a+ + (<i>„lV<i>_)a_], j = ±. 


We evaluate the inner products in (6.9) using the following two propositions. 

Proposition 6.1. 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


($+,i?2-Vx$-)^^(a) =0 , 

(4>_,1%2 • Vx4’+)^^(q) = 0 , 

2 ($+,1^2 • Vx4’+)^2^(.j 2) = +^l^#l l-^al 
2($_,i?2-Vx$-)^2^(n) = -^|A#| \^ 2 \ 


The constant, Aj e C, is generically non-zero; see Theorem 3.4 
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Proof. Let 32 = -^ 2 ^^ By (7.28)-(7.29) of 14 (see also 13 ) we have: 


;( 2 ) 


(6.14) 

(6.15) 

(6.16) 


($1,1^2 • Vx4>2)j;^^(f2) — 2 32 ) 

($2,1^2 ■ Vx<i)i)^ 2 ^(f^) = - Ajj 32 , 

($h, ^2 ■ Vx4>b)^2^(0) = 0i & = 1,2. 


Relations (6.10)-(6.13) follow from the expressions for <i)± in (6.6) and relations 

(|6lll)-(|6l6f~~^ □ 

Proposition 6.2. Assume that W{tc) is real-valued, odd and Ah— periodic. Let 
-dj = ($i,fL$i)j ^2 (Q). Then, 'djj € M and 

(6.17) = ($_,R^ci>+)^^(^), 

(6.18) (<i>+,VL$+)^^(f,) = ($_,R^ci>_)^^(^) =0 . 

Note that since $+(x) = e*^* ’'P+(x) and 4>_(x) = e®^*'’‘P_(x), relations (6.17) 
and (6.18) hold with $+ and $_ replaced, respectively, by P+(x) and P_(x). 

Proof. Equations ( |6.17 l-( | 6 )l 8 1 follow from the relations 

(6.19) ^ = - (‘i’2, , 

(6.20) (‘i>i,W^‘J>2)i^(n) = (<I’2,1R<1 >i)^^(^)=0. 


To prove (6.19) and (6.20), we begin by recalling that 4>2(x) = $i(—x), W is real¬ 
valued and 1T(—x) = — VE(x). Since W is real-valued, it is clear that dj S M. 
Furthermore, 

($ 2 , 1 E$ 2 )l 2 (o) = / ^' 2 (x)lT(x)$ 2 (x)dx = /" $i(-x)lE(x)$i(-x)dx 

' Jn Jn 

= j <i)i(x)lE(—x)$i(x)dx = — 111 ) . 

Jn 

This proves (6.17). To prove (6.18), observe that 

($ 2 , 1 T$i )^2 (Q) = [ $ 2 (x)lE(x)$i(x)dx = [ $i(-x)TT(x)$i(x)dx 

Jn Jn 

= I $i(x)lE(-x)$i(-x)dx = - / $i(x)lT(x)$i(-x)dx 


- -{^i,W^2)— ~ (‘^’2, bE4>i)^2^ 


in) 


This completes the proof of Proposition |6.2 


□ 


Propositions 6.1 and 6.2 imply that the orthogonality conditions (6.9) reduce to 
the following eigenvalue problem for a(^) = (a+(^),a_(C))^: 


( 6 . 21 ) 




Q — 0, o. ^ L 


Here, V denotes the ID Dirac operator: 

( 6 . 22 ) V = -i|Aij||i?-2|cr39c + i?iiK(C)o-i, and Aj x 7 ^ 0 . 

In Section 6.1 we prove that the eigenvalue problem ( 6.21[ ) has an exponentially 
localized eigenfunction a*(C) with corresponding (mid-gap) zero-energy eigenvalue 
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= 0. Moreover, this eigenvalue has multiplicity one. We impose the normal¬ 
ization: ||a*||i 2 (H) = 1. 

Fix = (0,a*). Then a* S is completely determined 

(up to normalization) and the solvability conditions (6.9) are satisfied. Therefore, 
the right hand side of (6.8) lies in the range of Hy — i?* : —> L^, and we may 

invert {Hy — i?*) obtaining 

(6.23) (x,C) + V'i^^(x,C) = + V'i^^(x,0, 

where 


R{E^) = {Hy - E,)-^ : 


P±H^ 


and P± is the L|^(n)— projection on to the orthogonal complement of the kernel 
of Hy — i?*, equal to span{$+,$_}. Here, denotes a particular solution, and 

V’i^^(x,C) = a+^(C)^+(x) +af!^(C)^’-(x) 

is a homogeneous solution. 

Note that by exploiting the degrees of freedom coming from the L^— kernel of 
Hy — i?*, we can continue the formal expansion to any order in S. Indeed, at 0{S^) 
for £ > 2, we have 

(6.24) (-Ax + y(x)-F;*)V’^^^(x,C) 

= ( 2(1^2-Vx) + 

+ (x, C; , ■ ■ ■, , 

-I- t), •) = •) V 0 e Ah, 

where, for the case i = 2, 

(6.25) 

G(2)(x,C;V'^°\V'^')) = ( 2(il2 • Vx) ac -^(Omx) )V'^iHx,0+|I^2p5^¥°Hx,0 . 


As before, (6.24) has a solution if and only if the right hand side is L^(n;dx)- 
orthogonal to the functions $j(x), j = ±. This solvability condition reduces to the 
inhomogeneous system: 


(6.26) 

(6.27) 


(C) = (C) + G 


where 


($_(.), gW (•, C; ,..., : E (^)))) 


L^m 


i^(n) 


Solvability of the non-homogeneous Dirac system (6.26) in L^(IR), is ensured by 
imposing L^(M)— orthogonality of the right hand side of (6.26) to Q!*(c)). This 
yields: 


(6.28) 




L2(R) 


Thus we obtain, at 0(6^), that = 4^^ + 4h \ where 4^^ i® ^ particular 




solution of (6.24) and 44 = a+HC)‘J^+(x) -I- (C)4>_ (x) is a homogeneous 

solution. 
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Summary: Given a zero-energy eigenstate of the Dirac operator, V (see 

Section 6.1), we can, to any polynomial order in 5, construct a formal solution of 
the eigenvalue problem i/; S ife||=K oi- 

6.1. Zero-energy eigenstate of the Dirac operator, V. 

Proposition 6.3. Let k(C) be a domain wall function (Definition \5.1}^ and assume, 
without loss of generality, that > 0. Then, 

(1) The Dirac operator, V, has a zero-energy eigenvalue, = 0, with expo¬ 
nentially localized solution given by: 


(6.29) 


a*(C) = 




— Ti-TT-T t^{s)dS 

g \><^\\ S2 \ -^0 ^ ‘ 


Eere, 7 S C zs any constant for which ||a*||j ;,2 = 1 . 


(2) The solution (6.29), a*, generates a leading order approximate (two-scale) 
edge state: 

-x) 


(6.30) 

(6.31) 


= a*,+ (( 5 f ^2 • x)$+(x) -I- a* _(( 5 f ^2 • x)$_(x) 
= 7 (-1-i) 






' K{s)ds 


(x, (5x) is propagating in the Oi direction with parallel quasimomentum 
A:|| = K • Oi, and is exponentially decaying, .^2 'X —> ± 00 , in the transverse 
direction. 

Proof of Proposition^^ The system (6.21) with energy E^^^ = 0 may be written 
as: 




0 1 

-1 0 


f a+ 
U_ 


and has solutions: 


^i(() = (^ 1 ) ei^ttii32i -^ 0 and / 32 (C) = f e. 


Since ilj > 0 and k(C) —t ±Noo as C —t ± 00 , with > 0, the solution / 32 (C) 
decays as C —t ± 00 . Thus we set Q!*(C) = 7 ,d 2 (C)j with constant 7 € C chosen so 
that ||q;*||^ 2 (r) = 1. This yields the expression for il'(°)(x, Jx) in (6.30)-(6.31), and 
completes the proof of Proposition |6.3[ □ 


Remark 6.4 (Topological Stability). The zero-energy eigenpair, (6.29), is “topo¬ 
logically stable” or “topologically protected” in the sense that it (and hence the 
bifurcation of edge states, which it seeds) persists for any localized perturbation of 
k(C). Such perturbations may be large but do not change the asymptotic behavior 
of k(C) as C —t ± 00 . 

7. Existence of edge states localized along an edge 
In this section we prove the existence of edge states for the eigenvalue problem: 
(7.1) = Evh , , where 

= — A -I- E(x) -I- (5k((5.^2 ■ x)IT(x) . 
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We make the following assumptions: 

(Al) 1^ is a honeycomb potential in the sense of Definition 2.4 and —A + V has 
a D irac point at (K, A*); see Definition 3.1 and the conclusions of Theorem 
3.4 In particular, the degenerate subspace of — A* has orthonormal 
basis of Floquet-Bloch modes {<l>i(x) , $ 2 (x)} and 

1 
i 






c(m)^ 


)(■ 


K + mk) 7 ^ 0; see (3.9). 


(A2) W is real-valued and A/j— periodic, odd and non-degenerate; ie. 
(W2) and (W3) of Section]^ hold. In particular, 

= 7 ^ 0 . 

(A3) k{5^2 • x) is a domain wall function in the sense of Definition 5.1 
The following spectral no-fold condition plays a central role. 


(Wl), 


Definition 7.1. [Spectral no-fold condition] Let Hy = —A -f I4(x), where D is a 
honeycomb potential in the sense of Definition | 2 . 4 | Further, let (K, A*) be a Dirac 
point for Hy in the sense of Definition | 3 . 1 [ in which we use the convention of labeling 
the dispersion maps by: k i-A- At,(k), where b € {6*, 6* -f 1} U {6 > 1 : 6 7^ 6*, 5 * -I- 1 } 
= { —1 : b —,+}■ 

To the Oi— edge, KOi, we associate the “.^2— slice at quasi-momentum K”, 
given by the union over all b G { — , -f} U {6 > 1 : 5 7^ —, -f} of the curves {(K -f 
Ai^2 , Bt,(K + ASl2} : |A| < i}. 

We say the band structure of Hy satisfies the spectral no-fold condition for the 
t)i— edge or, equivalently at the Dirac point and along the M2— slice, with constants 
Cl, C2, ao, and u G ( 0 , 1 ) if the following holds: 

There is a “modulus”, w(a), which is continuous, non-negative and increasing on 
0 < a < tto, satisfying 01(0) =0 and 

uj(a‘')/a—>■ 00 as a —)• 0, 


such that for all 0 < a < Rq: 

1 


(7.2) 

(7.3) 


a" < 


< 


6 7^±, |A| < 1/2 


E±(K + AM2)-E^ >cia;(a"), 

-£'b(K -f A.^2) ~ -S'* ^ C2 (1 -f \b\) . 


Our final assumption is 

(A4) — A-l-F satisfies the spectral no-fold condition at quasimomentum K along 
the M 2 — slice; see Definition 17.1 


Remark 7.2. (1) Conditions (7.2|-(7.3) ensure that, restricted to the quasi¬ 

momentum slice A I— >■ K -I- A .^2 £ ^h, the dispersion curves which touch 
at the Dirac point (K, A*) do not “fold over” and attain energies within 
Cl • a;(a'^) of A* for quasimomenta bounded away from K. 

(2) Dispersion curves of periodic Schrodinger operators on (Hill’s operators, 
H = —df Q{x), where Q{x -f 1) = Qix)) with “Dirac points” (see [9 12 ) 


always satisfy the natural ID analogue of the spectral no-fold condition 
with a;(a) = a. Dirac points occur at quasi-momentum k = ±7r and ODE 
arguments ensure that dispersion curves are monotone functions of k away 
from k = 0, ±7r. 
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(3) In Section we prove that H^y = —A -|- eV, where F is a honeycomb 
potential, satisfies the no-fold condition along the zigzag slice (Oi = Vi) 
with modulus a;(a) = a^, under the assumption that eVi^i > 0 and e is 
sufficiently small. 


We now state a key result of this paper, giving sufficient conditions for the 
existence of Oi— edge states of for Ui G A/j. 

Theorem 7.3. Consider the t)i— edge state eigenvalue problem, (O), where V (x), 
W(x) and k{Q satisfy assumptions (A1)-(A4). Then, there exist positive constants 
(Jo, Co and a branch of solutions of (O), 

|J| G (0,(5o) ^ G (£;* -coSo, E^ + co Jo) x 

such that the following holds: 

(1) is well-approximated by a slow modulation of a linear combination of de¬ 
generate Floquet-Bloch modes ih-i- and $_ ( (6.6 )), which is decaying trans¬ 
verse to the edge, ZOi; 


(7.4) 

(7.5) 


¥{■) - K+(Jj^2-)®+(-)+a*,-(^-S^2-)^-(-)] 




< J5 


E^ = E„ + + o((j2), 

where E^"^^ is obtained directly from (6.28), (6.27) and (6.25). The implied 
constant in (7.4) depends on V, W and k, but is independent of 5. 

(2) The amplitude vector, a*(C) = (a*_-|-(^),(<C)), is an L^(IR^)— nor¬ 
malized, topologically protected zero-energy eigenstate of the Dirac system 
(6.22).- Pa* = 0 (see Proposition \6.Sl\ l. 

Perturbation theory for fc|| near K • t)i can be used to show the persistence of 
edge states for parallel quasi-momenta near K • tJi. 


Corollary 7.4. Fix V, W, k and 6 as in Theorem \7.3\ Then there exists pq S 
such that for all fey satisfying |fc|| — K • t)i| < ryo, there exists an (S)— eigen¬ 
function with eigenvalue /i(J, fey) = E^ -\- pf (fcy — K • Oi ) -fO(|fcy-K.0i|2), 
where E^ is given in (7.5), and is a constant, which is independent o/fcy. 


Zigzag edge states for /cy in a neighborhood of /cy = K • Oi = 27r/3 (Oi = Vi) and, 
by symmetry, in a neighborhood of fcy = 47r/3 are indicated in Figure 

7.1. Corrector equation. We seek a solution of the eigenvalue problem ( |7.1[ ), 
in the form 


(7.6) 4-^ = i/((°)(x, SSi 2 ■ x) -k ■ x) -k 6 rj\x.), 

(7.7) E^ = E,^ + d^p^. 


Here, and are given by their respective multiple scale expressions ( |6.7[ ) 
and ( |6.23[ ): 

V'^°^(x,(Jil 2 • x) = a*,_|_((Ji ?2 • x)$+(x) -k a*,_((Ji ?2 • x)$_(x). 


i/i(^)(x, Jx) = (x, Jj ^2 ■ x), 

and {p^,ri^{x)) is the corrector, to be constructed. We may assume throughout 
that (J > 0. 
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Remark 7.5. We shall make frequent use of the regularity of <f>+(x), $-(x) and 
a*(^) = (a*^+(C),a*,-(C))^- In particular, V G C''^(M^/A) and elli ptic regularity 
theory imply that is (^““(M^/A), and by Proposition 6.3 a*(C) and its 


derivatives with respect to C are all exponentially decaying as |C| —t oo. 

The following proposition lists useful bounds on and . 

Proposition 7.6 (77 ^||=k.oi(^x) bounds on (x, 5.^2 • x) and V'p^^(x, 5.^2 ■ x)). 
For all s = 1,2,..., there exists Jq > 0, such that if 0 < |(5| < (5o, then the leading 
order expansion terms SSi 2 ■ x) and ' x) displayed in (6.7) and 

(6.23) satisfy the bounds: 


i/;(°)(x, SSi2-^) 


It follows that II //2 


Hi 




(^=(5.^2-X 


V'^^^(x, SSi 2 ■ x) 

9.9cV'W(x,C) ^ 

C=d.H2-X 


Hi 




= C>((5i/2). 


The proof of Proposition 7.6 follows the approach taken in the proof of Lemma 
6.1 in Appendix G of 12 . We omit the details but make two key technical remarks, 
that facilitate this proof. 

Bound on ||$h|li/».- Recall $ 6 (x;k) = e**"'^p{,(x;k), where 11 ^* 61 ^ 2 ( 0 ) = lbb||L 2 (n) = 

1 . Now p 6 (x;k) satisfies —Apt = 2ik ■ \7pi, — \k\‘^pb — Vpb + Eb{k)pb, where V 
is bounded and smooth. By 2D Weyl asymptotics |£’ 6 (k)| « (1 + | 6 |), 5 » 1 
and therefore we have ||Ap 6 ||h»-i < Cs(l + b) IbbU/i'’- Hence, by elliptic theory 
Ibbllffo+i < Gs(l + b) IbblU" and induction on s > 0 yields \\pb\\H‘ < Gs(l + 6 )® . 

Rapid decay o/($ 6 (-; K),/(•))i 2 (n) • Using ij(o)$ 6 (x;K) = L; 6 (K)$ 6 (x; K), for suf¬ 
ficiently smooth / we have: ($ 6 (-;K),/(-))i 2 (n) = (U 6 (K))-i'^($ 6 (-;K),[i 7 (°)]i'^/(-))^ 2 (n)- 
Hence, for any M > 0, I ($ 6 (-;K),/(-))i 2 (f 2 ) I < CM{^ + b)~^. 

It remai ns to construct and bound the cor recto r (/i,77(x)). Substitution of the 
expansion ( |7.6[ ) into the eigenvalue problem (7.1), yields an equation for r)(x) G 
'Ufc||=K oi (U)? depends on p, and the small parameter <5: 

(—Ax -I- H(x) — A*) ? 7 (x) + 5k{S^2 ■ x)lT(x)? 7 (x) — S'^p j)(x) 


= S 
(7.8) 


(" 2.^2 • Vx b — K(( 5 .ft 2 • x)lT(x)')'i/'b^(x, C)! + 5 mU*‘°^(x, ( 5.^2 • x) 

V / Ic=<5A2'X 


-I-(5|1^2|^9cV'^“’(x,C) 


C = 5j?2-x 


-I- (51^2 • x) -I- ( 5 ^| 1 ^ 2 |^ dc V'p^^(x, C) 


(^ = S^2 'X 


7.3 


To prove Theorem 
G ^K oi satisfying the bound 


we shall prove that (7.8) has a solution {p{6),p^), with 

< C'<5U^ . 
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7.2. Decomposition of corrector, 77 , into near and far energy components. 

Introduce the abbreviated notation, for |A| < 1/2: 


'E,iK + X^2) K^{K,K + 1}, 


(7.9) 
and 

(7.10) 


Et{X) = {E_i\) 

E+iX) 


b = K, 

b = bi, + 1, 


‘h 6 (x; K -|- X^2) b-i, ^ {&*, 6 * + 1 }, 

$b(x;A) = <($_(x;A) 6 = 6 *, 

$+(x;A) 6 = 6 * + l. 


Define /h(A) = ($ 6 (-, A),/(•))l 2 . By Theorem 

)t||=Koi 

has the representation 


4.2 


any 77 e 77^ k.oi(S) 


(7.11) 




$f,(x; A) ? 7 b(A) dX . 


b>i •>'DI<i/2 

Our strategy is to next derive a system of equations governing { 77 b(A)}b>i, which 
is formally equivalent to system (7.8 1 . We then prove this system has a solution, 
which is used to construct ? 7 (x). 

Take the inner product of (7.8) with $h(x; A), for 6 > 1, to obtain 
6 > 1 : (EbiX) - E,)?jb{X) 

(7.12) +(5($b(-;A),N(6l^2-)m-)^(-))L= ^ 

fc||=KDl 

= 6Fb[t^,6]{X) + 6^ fi?jb{X) , |A| < 1/2. 

Here, (5](A), 6 > 1, is given by: 

(7.13) Fb[l^,S]{X)^F^’\X) +t^F^’\X) + 6^F^’\X) + F^’\X) + 6F^’\X), 
where 

F/’"(A) = /<E>6(x, A), (21^2 • Vx dc - <SSi2 ■ x)lT(x))7/;/i)(x, C) 


C = 5 .^ 23 c/ 7-2 


F/’^(A) = (4),(x,A),V>'"J(x,5l^2-x) 
(7.14) 

Ei’HX) = (-I>i,(x, A),V)p^'(x, 5 l ^2 • x) 


F/-^(A) = ^-l>,(x,A), |i^ 2 |' 9c^'°Hx,C) 
Eb’\X) = ($i{x,X), |J^ 2 |^ C) 


C = 5.^2-x/ j-2 

^ku =K-Oi 


C= 5 .^ 2 -x/ 7-2 

fc II =K • o 1 


Recall the spectral no-fold condition ensuring that 5/oj{5'^) —>■ 0 as 6 —?> 0, where 
p > 0. We next decompose 77 (x) into its components with energies “near” and “far” 
from the Dirac point: 

77 (x) = 77„ear(x) -|- 77far(x), where 
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(7.15) r?near(x) = ^ 

b=± 

(7.16) ??far(x) = ^ 


/|A|<l/2 


$6(x; A) 77b 

,near (A) d\ 


b>l ■^|A|<l/2 
^±,near(A) = x(|A|<5‘') 77±(A), 


$b(x; A) 776,tar(A) dA, and 


77t,jar(A) = X +(5,,_)(5^ < |A| < -j 77t,(A), 6 > 1; 

<5b.+ and 5h,- are Kronecker delta symbols. 

We rewrite system (7.121 as two coupled subsystems: a pair of equations, which 
governs the near energy components: 

(if_|_(A) — if*) 77 +„ear(A) 

(7.17) + 5x[ |A| < <5“^) ($+(•, A), «(di?2-)W^(-) [77near(’) + dfar (’)] )l2 ^ ^ ( 2 ) 

= Sx{ |A| < d'^)F+[/X,<5](A) + dV d+.near(A), 

{E_{X) - i;*)77_,„ear(A) 

(7.18) +dx(|A| < S'') ($_(•, A), k(( 51^2-)W^(-) [?7near(-) + dfar(-)])L 2 

fc|l 

= <5X(|A| < S'')F_[fi,S]{X) + dV d-.near(A), 
coupled to an infinite system governing the far energy components: 

(ifb(A) — if*) 77byar(A) + Ay^1/2 ^ IA| > (Sb^- + 6'')x 

(7.19) ($b(-, A), 7t(A.ft2’)^(') [dnear(’) + dfar(’)])Lj (S) 

= Ax(l/2> |A| > (Ab,_+Ab,+)A'^)Fb[/i,A](A) + S 17 b,far (A), b 


,(S) 


> 1 . 


We now systematically manipulate (7.17)-(7.19) into the form of a band-limited 
Dirac syste m; see Proposition |7.12| This latter equation is then solved in Propo¬ 
sition 7.15 Since all steps are reversible, this yields a solution ( 77 ^, {7?b(A)}b>i) of 
(7.17)-(7.19). Finally, r]^ G solution of corrector equation (7.8), 

is reconstructed from the amplitudes {77b(A)};,>i using (7.11). 

7.3. Construction of 77far = 77far[77near, A] and derivation of a closed system 
for 77near- We solve (7.19) for 77far as a functional of 77„ear, and the parameters /i 
and A. We then study the closed equation for 77near obtained by substitution of 77far 
into (7.17) and (7.18). 

It is in the construction of this map that we use assumption (A4), the spectral no¬ 
fold condition along the ^ 2 — slice; Definition \7.1\ We apply it in the form: There 
exists a modulus, a;(a), and positive constants v, ci and C 2 , depending on V, such 
that for all S 0 and sufficiently small: 

(7.20) ^ 


S'' < 


< 


(7.21) 


6 7^ ± : |A| < 1/2 


E±{X)-E, > ciw(A"), 

EbiX)-E, > C2 (1 + |6|) . 























38 


C. L. FEFFERMAN, J. P. LEE-THORP, AND M. I. WEINSTEIN 


The far energy system (7.19) may be written as a fixed point system for lyfar = 

{ ^b,far ('^) }&> 1 • 

(7.22) 

: ^near 5 — ^b,far ) ^ ^ I5 

where the mapping £’& is given by 

?Mar(A) Sx{l/2>\M>iS,,-+6b,+)6’') 




Eb{X) - E, Eb{X) - E, 

- {M; A), ^iSSi2-)Wi-) m + m)Ll _ + Fbifi, ^](A) ) , 


and 

'/'(x) = / x( |A| > {Sb- + <5&.+)(5'') (^h(A)$ 6 (x; A) dX 

= y] / ^&yar(A)$b(x; A) dA . 

5>1 

Equivalently, 

(7.23) ^[^far5 ^neari/^;— ^far ■ 

For fixed S and band-limited 77near: 

(7.24) ^±.near(A) = X ( | A| < (5'") ^±,near (A), 

we seek a solution {^t,,far(A)}h>i, supported at energies bounded away from if*: 

(7.25) i)'t, far(A) = X(|A| > (dfe- +db,-|-)(5'^)^6,far(A), 6 > 1. 

Introduce the Banach spaces of functions limited to “far” and “near” energy regimes: 

iLar.^ys) = {/ e LI||.k.„,(S) : A(A) satisfies , 


= {/ e tI„=k.„,(S) : h{X) satisfies ([g} , 


Near- and far- energy Sobolev spaces iif^^(S) and are analogously defined. 

The corresponding open balls of radius p are given by: 


-2 

''near,(5 


< P 

-^fc||=KDi 


Tinear.d^ (p) — ^ ^ 7/ 

Bfar,S’'{p) = |/ G < p| • 

Using (A4) that = —A + V satisfies the no-fold condition for the t>i— edge, 
we deduce: 

Proposition 7.7. (1) For any fixed M > 0,i? > 0, there exists a positive num¬ 

ber, So < 1, such that for all 0 < S < (5o» equation (7.23), or equivalently, 


the system (7.22), has a unique solution 


{Pnear-; P: S') G linear,(7?) ^ {|p| ^ iff} X {0 < d < do} 

' ^ dfar(‘7 bnear? P: d) = E Pfar G Efar^S^ iPs) ^ P5 ~ E) 


ds 

a;(d^ 
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(2) The mapping (rinear, pl, S) 77far(-;’7near, <5) e ■f^fc||=K Oi Lipschitz ir 
(^near,M) with: 

II VfarblJl,Tl,5]-'qfar[lp2,t^2,S] \\^2 

6 


< c 


UJ(S‘') 

(7.26) II TlfarlVnear', ||//2 


( ll^/’l - V’2||//2 +\tii-H2\), 

\ fe||=KDl / 

<5 „ 


< C" 


l inear f72 + /c,.\ 


T/ie constants C and C" depend only on M,R and v. 

(3) The mapping (r?„ear, <5) H> PfariPnear, hi',S\ satisfies: 

(7.27) Pfar[p near-) 5](x) = [Ap 

near ](x; S) + nB{yL-, 5) + C'(x; <5). 
For Pnear G .Bnear(-R) WC haVC: 

6 


mp near] (‘; ^1,(5) - [Ap 

near ](-,M2,^)||//2 
fc 

5 






Imi - M 2 I 


II near ](•;Ai,^)li Hi 


— / r7/\ K?near ff2 

fc|l=Kt.i Uj{0'^) "-"fc||=KDi 


l|i?(-;<5)ll- 


< 


and ||C'(-; 5)11^2 


< 


(4) We may extend ?7far[’; ^nearj <5] to be defined on the half-open interval S G 
[0,(5o) by defining 77far[?7near, <5 = 0] = 0. Then, by ( |7.26[ ) 77far[?7near, M, <5] 
is continuous at S = 0. 


Remark 7.8. [Remarks on the proof of Proposition 7.7 The proof follows that of 
Corollary 6.4 in [^, with changes that we now discuss. 

(a) The fixed point equation (7.22) for ■i/'far is of the form: 


(7.28) 


, <5 x(|A| > ((56.-+ 56.+)<5‘^) 

m.r - QsP!.r + EbiX)-E^ 

-{MGX),R{SSi2-)Wi-) VneA-))Ll 


F,[pi,d]iX)], 


where Qs is bounded and linear on ff^n^K-oi defined by: 


W]sW ^ -S (ci>6(-,A),^(,5jl2-)5y(-) 


(7.29) 


S^pL 


Et,{X) - E, 

MX) 


EbiX) - E, ■ 

To construct the mapping (pnear, P-, d) 1 —77far [??near > P, and obtain the con¬ 
clusions of Proposition [ 7 ?^ it is convenient to solve ( 7.28| ) via the contraction 
mapping principle. Thus we need to bound the operator norm of Qs and 
we find from ( |7.29[ ) that Qs maps Tnear.a-' Mear,S’' norm bounded 
by constant x e((5), where 


(7.30) 


e(5) = sup sup 


l<5| 


*'=± I'Sp^lAl^i \^b[A) --H*\ b>l, b^± 


+ sup ( 1 -b |6| ) sup 


1^1 


o<iAr<ii^^(W-^’^r 
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(7.31) 


The spectral no-fold condition hypothesis (7.20)-(7.21) implies that 

1^1 


e(h) 


< 






w(5-) CiiV) C2{V) ’ 

which tends to zero as <5 tends to zero. Hence, the contraction mapping 
principle can be applied on the ball B{ai-^s‘'{ps), ps = 0{S^/oj{5'')). 

(b) We note that although E is a two-dimensional region, since E is unbounded 
in only one direction, estimates on (E) have the same scaling behavior 
in the parameter 5 as in the ID study |l^ . 

7.4. Analysis of the closed system for lynear- Substitution of ??far[?7near, P: <5] into 


the system (7.17)-(7.18) yields a closed system for (?7near,M)i which depends on the 


parameter 5 S [0, 5o). In this section we show, by careful rescaling and expansion 
of terms, that the equation for pneai may be rewritten as a Dirac-type system. We 


then solve this system in Section 7.5 Recall the abbreviated notation: i?&(A) and 


$f,(x;A), introduced in (7.9)-(7.10). 

Since both the spectral support of ??near (parametrized by K -|- with |A| < 
5''), and size of the domain wall perturbation, 0{5), tend to zero as 5 —>■ 0, it 
is natural to scale in such a way as to obtain an order one limit. We begin by 
introducing a scaling of the quasi-momentum parameter. A, and 77±,near (Oj an 
expression for 77±,near(A) as a standard Fourier transform on K: 

(7.32) ^it.near (0 = ^ib,near(A), where —. 


By Proposition 4.5 (A)—if* = ± |Aj| |.ft 2 | i^^+if 2 .± (<^0 (J^)^, where |. E2.±( i5p| < 

1, for all s ee (4.12). Substitution of this expansion and the rescaling (7.32h into 
(7.17)-(7.18), and then canceling a factor of 5 yields: 


(7.33) 


+ I AbI IA 2 U 57+.near(e) + Xd?! < {<!>+(■, 5^ , ^(<5^2 • ) VF(- )hnea. (■) ) i 

= XdCI < 5](A) + 5/4 h+,nea.(C) - 5^2,+(50C"h+,near(0 

-xd^l < 5'""^) {4>-H(-,5C),N(5A2-)IF(-)4?far[l7near,/i,5](-))i2 


- |Ab| IA 2 I e 57-,near(C) +X(I^I < 5"“') (4> _ (•, 5^, ^(5 A 2 •) IF(-)hnear (•)) L? 

'=11 

(7.34) = xdCI < 5"-l)i?_[M,5](A) + 5^ h-,near(0 - 5i?2,-(50C"h-,near(0 

-Xd^l < 5"-'){$-(-,50,N(5A2-)IF(-)hfar[l7near,/4,51(.))i2 


We next extract the dominant behavior, for 5 small, of the inner products involving 
i?near by first expanding r/near in terms of its spectral components near energy if* = 


if±(A = 0) plus a correction. To this end we apply Proposition 4.5 to expand 
p± (x. A) for \ = 5^ small: 


p±(x,A) = P±(x) -f (/j±(x,5^), P±(x) = ^ |^|^Pi(x)±P 2 (x) 


where 


(7.35) 


|7'±(x,( 50I< sup |(/?±(x,a;)| < d"", |C| < ^ 

xGS, 


Thus, using (7.15) and that 4)±(x; A) = x A) (see (4.11)), we obtain 
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r (^) — / <E>+(x,A)j7 + ,near (A)dA + / <1 ?_(x,A) 77 —,near iX)dX 

J\x\<5^ JwKS’^ 

-L 


'W<S' 

/ 

f\\\<5- 

I 

f\X\<5' 


A)77_,near ( ^ ) dA 


(7.36) = Se'^' 


iK x zAAs'X / 

e e p+(x,A)77+,: 

I iKx i. 

/ e e 

V|Ai<5^ 

(x) [ + ^ 6 *^ "‘P+(x, ' x) 

^-(x) [ e“*^^^'"‘77^_.near(C)rfC + *P-(X, (5i^2 • X) 

d|5|<5—l 

^’+(x) P+ ,near (d.i ^2 ' x) ~t“ d^_ (x) 7^—,near(d.i ^2 ' x) ~t~ ^ ) Pb(x, d .^^2 ' x) 


r iKx j-j 

= oe P- 


+ de*"P-(x) 


(7.37) where p±(x, C)= f e^‘’ifi±{x,S^)rj±^ne3.T{0dC 


6=± 


We now expand the inner product in (7.33); the corresponding term in (7.34) is 
treated similarly. Substituting (7.36) into the inner product in (7.33) yields (using 
^’±(x;k) = e*'"’"p±(x;k)) 


(^’+(-,^C),K(l5-)W^(-)»7near(-))L2 

(7.38) = ($+(•, K + <5ei^2),«((5-)W^(-)r/„ear(-))L== 

ft II =K-tii 

(7.39) = W^(-) ^+,near(<5i^2-))i2(s) 

(7.40) +6{e^^^^^-p+{;dO,P-{-) H^(-) <S^2-) 77-,near(^i^2-))i2(s) 

(7.41) +5^ p+i;60,W{-) KiSSi2-) Pb{;5^2-))^,^^^ ■ 


The inner product terms in (7.39)-(7.41) are each of the form: 


(7.42) Q{6;^) = 6 f e (5^)r(x, < 5.^2 • x)dx, where 

JT. 

(IPl) g{x,y) is a smooth function of {x,y) G x M and 

(IP2) x hG r(x,() is A/j— periodic and with values in L^(IR^), i.e. 

(7.43) r(x + v, C) = r(x, C), for all V G A/i, 

(7-44) H H / ll^xr(x,C)lli2(Rj dx < 00. 


We denote this Hilbert space of functions by with norm-squared, || • 11 ^ 2 , 
given in (7.44). It is easy to check that conditions (7.43)-(7.44) are satisfied 
for the cases P = r(C) =J<C)??±,near(C) and P = P(x,C) = «:(C)p±(x,C), 
where p± is defined in (7.37). 

To expand expressions of the form G{S,^), we use: 


Lemma 7.9. Let g(x,y) andP(x, (^) satisfy conditions (IPl) and (IP2), respec¬ 
tively. Denote by P(x, w) the Fourier transform ofT{x,f) with respect to the f— 
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variable, given by 

(7.45) f(x,a;) = lim ^ f 

Atoo Ztt J\(^\<n 

where the limit is taken in x dxda;). Then, 


(7.46) 


nGZ'^^ 


n 




with equality holding in Lfo(.([“^max, ■fmaa:]; rfC)> any fixed > 0. 

We adapt the proof in (Lemma 6.5) for the ID setting. We require the following 
variant of the Poisson summation formula in 


Theorem 7.10. LetT{x,() satisfy (IP2). Denote byT{x,uj) the Fourier transform 
o/r(x, C) with respect to the variable, C; see (7.45). Fix an arbitrary i/max > 0, 
and introduce the parameterization of the cylinder Ti: 'k. — TiXti T 2 X) 2 , 0 < ri < 
1 , T2 € K. Then, 


E 6 til + r 2 t) 2 , T 2 + n) = 2-71 E (riPi + r 2 p 2 ,27rn -|- y) 

n^'L n£Zi 


in ([ 0 , 1]^ X [-yinax, ymax]; ^^1^x2 ’ dy). 


. Since the 

proof is very similar, we omit it. We also require 


The ID analogue of Theorem 7.10| was proved in Appendix A of 12 


Lemma 7.11. LetF{:>c,y) andFNfK,y), N = l,2,..., belong to L'^{T,x[-y^,,,^,y,nax]]dTcdy). 
Assume that 


II^A --F|lL2(Sx[-y„,,,,y_];dxdy) ^0^ SS N ^ OO . 

Let G G L^([—j/max, 2/max]; dy)- Then, in the L'^{T,; dx.) sense, we have: 

^2/max Pyrciax 

lim / FAf(x,?/)G(y) dy= lim Fjv{x,y)G{y) dy 

N^OC ' I AT . . . 


'-2/rr 


' Umax 
/* J/max 


N—^oo 


F(x, 2 /)G(y) dy. 


Proof of Lemma 7.11\ Square the difference, apply Cauchy-Schwarz and then inte¬ 
grate dx over E. □ 


Proof of Lemma [ 7. Recall the parameterization of the cylinder, E: 

X G E : X = TiPi - 1 - T2X12, 0 < Ti < 1, T2 G K , .fti • x = 2 -kti, R2 ■ ^ = 27 rT 2 , 

dxi dx2 = jOi A 02 ] dri dT2 = |n| dri dT2. 


Using that y(x, 5^) = g{TiXii F T 2 V 2 , S^) and r(x, 6^2 ■ x) = r(r iPi 4 - T 2 O 2 , ‘ 2 .'k 5 t 2 ) 
are both appropriately 1-periodic, we expand defined in (|7.42 ). By Lemma 
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[ml 

G{5-,0 

= 5 |Q| 

Jo 


= 5 |n| 

r dr, 

Jo 


= 5 |n| 

r dri 
Jo 


= 5 |n| 

rdri 

Jo 


n= — N '' 
N 


By Theorem 7.10 with T = r(rit)i + Tt) 2 , 2 ttSt2 ) and y = we have 
^ ^ TX> 2 ,2^5X2) = J I] (nOi + TX> 2 . ^ + C) . 

n^Tj n^'L 

with equality holding in ([0,1]^ x [—^max) Cax]; dTidT 2 ■ d^). Again using Lemma 
|7.11| we may interchange the sum and integral to obtain 

G{S;0 = ^mJ dn J 5(tiOi + r2t)2,i5C) +r2t)2, ^+C) 

(7.47) 


nGh 


= Y, [ (x, ” + e) 5(x, SOd^ . 

riGZ"^^ 


This completes the proof of Lemma 7.9 


□ 


We next apply Lemma 7.9 to each of the inner products (7.39)-(7.41). 
Expansion of inner product (|7.39 |; 

Let g(x,^^) =p+(x,i5^)P+(x)lT(x) and r(x,C) = K(C)7?+,near(C)- By Lemma 
S{E^^^^-p+i-,S0,PAW{-HS^2-)V+,ne.ASSl2-))i^.^^^ 


7.9 


= E 

n^Z *■ 




[k?7+ 

,near] ^ j + C) 7 ’+(x, (5C)P+(x)W(x)dx. 


Since p±(x, A = S^) = T’±(x) + (p±{x,6^), where (p±{x,6^) satisfies the bound 


(7.35), we have 

6 {E^^^^-p+{-, 60 ,P+{-)W{-)<S^ 2 -)V+,ne.rm 2 -))^ 2 ^^^ 

= ^+,near) + d+i^i V+,near), 


where 

(7.48) 

U(^’^+.near) = ^ ^ d'i 


nearJ — / . (^[^^+,near 

nGZ 






p+(x)rw(x)dx, 


(^;??+.near) = Xl-^cN+,near] +C) f S^)P+{x)Wix)dx. 

riGZ 


From Proposition | 6 . 2| and Assumption (W3) we have 

(7.49) f |P+(x)|^ lF(x)dx = 0 and f P+(x)P_(x)lF(x)c?x = 7 ^ 0. 

Jn Jn 
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Therefore, the n = 0 term in the summation of 77 +,near) in (7.48) is zero and 

we may write: 


4 


(C;»7+,near) = -^cN+.near] + C) ^ |P+(x) W^(x)dx. 


|n|>l 


Expansion of the inner product (fflot.- 

Similarly, with g{:ic,Sf) = p+(x,5^)P_(x)W^(x) and r(x,^) = k{C)p- 

,near (C), we 

have 

= 7^ 77—,near) T77_ near) 5 


where (noting, by (7.49), that the n = 0 contribution is nonzero) 

'^4-(Ci ,near) = 'du/v?7_|_ near(0 4“ A .near ) j where 

/+(C;?7-,near) = X! (l'^^) / A-(x)n"(x)dx, and 

/RC;?7-,near) = AN-,near] +^) j 6*"^" (x, A)-P-(x)!^ (x)dx. 


nel 


Expansion of inner products (7.41); 


Consider the h = + term in (7.41). Let g{-x.,5f) = p+(x, A)W^(x) and r(x,(^) = 
k(C)p+(x, (^). By Lemma 7.9 and the expansion of p+(x, 5f) about P+(x) in (7.35) 
we have: 

5 (e*«^-^=>+(-, A), W^(-)k(<5A0P+(':<5A-))l2(s) = ^+(?;^+.near) +/+(?; ?7+.near), 
where 


I+(?;^+.near) = [ NN+] (x, T + ^] 

nGZ“'^ 

/+(?;'?+,near) = ^ [ AN+] (x, T + ?^ 

e*”-‘^=“'^P+(x)lF(x)dx, 

g*nii2-x^_^(x,J^)^(x)dx. 

For the b = — term in (7.41) we have 


J (e*«^^-p+(-, A), VF(-)«(^A-)P-(-, <5A-))i.(s) 

= /+(?;r?-,near) +^+(?;i7- 

where 


/+(?; 17 +,near) = X! / NN-] (x, T + ?^ 

/+(?; 17+,near) = [ NN-] (x, ^ 

nGZ'^^ 

e“-^=“xp^(x)lF(x)dx, 

e“-^^'X^+(x,A)W^(x)dx. 


Assembling the above expansions, we find that the full inner product, (7.38), may 
be expressed as: 

8 

(7.50) ($+(•, A), K(<5A-)W^(-)»?near(-))L2(E) = A Anear,-(?) + ’ 

1=1 
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A similar calculation yields: 

8 

(7.51) ($_(•,(50,K((5l%2-)W^(-)?7near(-))L2(E) = '^t^near,+ iO+'^lU^lVnee.r), 

i=i 


where the terms 7i(0^near) are defined analogously to /^(Oi?near)- We now sub¬ 
stitute our results (7.50)-(7.51), (7.13)-(7.14) and (7.27) into ( |7.33 1 -( 7.34 1 to obtain 
the following: 


Proposition 7.12. Let = (77+,„ear(0,^-.near(O)^- Equations ( |7.33| )-( [7)34| , 
the elosed system governing the near energy components, ? 7 near, of the corrector, rj, 
is of the form: 


7.331-( 7.34 


(7.52) (V^ + C^ip) - Sfa) m = 6) 6). 


Here, denotes the band-limited Dirac operator defined by 


(7.53) = |A#| | 1 ^ 2 | CT3 ^ m + x(IO < txi nPiO- 


The linear operator, C^{fT), acting on f3, and the source terms Ai{^;6) andJf{^;S) 
are defined by: 


(7.54) 13(f) = X {\f\ < S'' '‘^here 


C{(g)m 

Ei(p)m 


Sf 


2 f E2,+ {Sf) 77 +,, 


rif) 

E2,-{5f) r7-,near(C) 
{<i>+(;Sf),KiSii2-)W(-)[Ag, 
{‘t>.(;5f),n(5Si2-)Wi-)[Av, 


Eiipfm = E 

Lear](-; M, ^))l2 „ 

fe||=K 

iear](-; P,'5))2^2 


I+(£,\g±,ne^r(f)) 
IL (C;?7±.near(0) 


3 

(7.55) A1(C; 5) = X (ICI < 5'^”^) ^^A4j(f; (5), where (inner products over 7/fe||=K.oi) 

i=i 


Mi(e;5) 


M2{f-,5) 




M2{f\S) = 5 


({^+(;5f),K(5^2-)W(-)B(--5))\ 

V{4._(.,,50.n(<51^2-)W(.)B(-;5)V ’ 


//$+(•, (50, 


4 

(7.56) 5) = X (ICI S’'~^)'^^Mj(f\5), where (inner products over 
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m^;S) = 


\ 


|^3>-i-(x, 50, ( 2 i ^2 • Vx 9c - ii{S ^2 ■ x)W(x)) ■ x)^ 

^$_(x,50, ( 2-^2 • Vx 9c — k{S ^2 ■ x)M^(x)) i/)p^^(x, 5,^2 • x)^ 

^^+(x,50, l^ 2 p 9^V'^°’(x,C) 

'$_(x,50,|^2p 9^V°nx,C) 

' ^+(x,50, l^al^ d^ipi^\x,C) 

$_(x,50,l-*^2p 9^1^^^^(x,C) 



Q- (f.x\ = _( ('*’+(x-^ 0 ,'^( 5^2 •x)W"(x)C(x; 5 ))\ 
l^(<i,_(x,50,N(5i^2-x)W^(x)C7(x;5))J ' 

We conclude this section with the assertion that from an appropriate solution 
(^f3^, fi{6)^ of the band-limited Dirac system (7.52) one can construct a bound 
state of the Schrodineer eigenvalue problem (|7.1|). We say f S jf 


Proposition 7.13. Suppose, for 0 < S < Sq, the band-limited Dirac system (7.52) 


has a solution y/3°(f,) , p{S) j , (3° = (/3+,/3_)“^, where suppj5° C {|0 0 satis¬ 

fying: 


< 




S 0 < S < 6o {to be verified in Proposition 1 7.7 5[) , 


p{S) bounded and p,{S) — p,o —>■ 0 as d —>■ 0 {to be verified in Proposition 1 7.7 9| ) . 
Define 

(7.57) flie.r,+ {0=Pii0, fne.r,-{0=P-{a 

and construct -I- as follows: 

(7.58) 






'/near,6 


$b(x; A)dA, 


^L,b(^) = ’7far.b[i?near,7‘, b > 1] (scc Proposition^^, 


9fa 


r(x) = H 


»7fa 


_ r.b(A)$b(x;A)dA+^ (A)$f,(x;A)(iA. 

;,=± •^'5‘'<|A|<1/2 f,^±"'DI<l/2 

?7^(x) = f?iar(x) + ??L(x), = E^, + d^p{d), Q<d <do. 

Then, for all 0 < |5| < 5o, the following holds: 

(a) i?^(x)ei7|„^K.o,(S). 


(7.59) 


(b) {ri^,p{d)) solves the corrector equation (7.8). 

(c) Theorem 7.3 holds. The pair {"i'^,E^), defined by (see also (17;§-(173; 

(x) = (x, X) + (x, X) + dp^ (x) , X = 51^2 • X, 

E^ = i?* -|- + o{d'^), 


is a solution of the eigenvalue problem (El with corrector estimates as¬ 
serted in the statement of Theorem\7.3\ 
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To prove Proposition |7.13| we use the following lemma. 

Lemma 7.14. There exists a (5o > 0 such that, for all 0 < S < Sq, the following 
holds: Assume /3 G L^(M) and let be defined by (7.57l-(7.58|. Then, 

< 5^/^ I 


l^near||//2 


lL2(R) • 


The proof of Lemma 7.14 parallels that of Lemma 6.9 in 12 , and is not reproduced 
here. 


Proof of Proposition \7.1t^ From ^ we construct such that: ||? 7 near ||//2 


< 


(5^/^ ||,^||j;, 2 (]i{) (Lemma 7.14). Next, part 2 of Proposition 7.7 (7.26), gives a bound 
on 77far: Ijiyfarliynear; (5] II//2 ^ C ( u]{S'') near 11^2 + These 


<5 2 

u(5-) 


two bounds give the desired Hf ^k.oi(^) bound on r]°. Note that all steps in our 


derivation of the band-limited Dirac system (7.52) are reversible, in particular our 
application of the Poisson summation formula in Therefore, given 

by ( [L^ is an Hi =k.b eigenpair of (|7.1[) . □ 


We focus then on constructing and estimating the solution of the band-limited 


Dirac system (7.52). 


7.5. Analysis of the band-limited Dirac system. The formal d 0 limit of 
the band-limited operator T>^, displayed in (7.53), is a ID Dirac operator V defined 


via: 


(7.60) 


vm = iabi ii^2i ^m)- 


Our goal is to solve the system (7.52). We therefore rewrite the linear operator 


in equation (7.52) as a perturbation of V (7.60), and seek /3 as a solution to: 
(7.61) vm +[&-V + C\fi) - Sfi) m = TM{i-, 5) + 5). 


We next solve (7.61) using a Lyapunov-Schmidt reduction strategy. By Proposition 
6.3 the null space of T) is spanned by a*(^), the Fourier transform of the zero energy 


eigenstate ( |6.29[ ). Since Qf*(C) is Schwartz class, so too is a*(^) and 3*(^) G iJ® 

), introduce the orthogonal projection operators, 

-^ 11 / = /)l2(h) S*, andPi/ = (/-P||)/. 


for any s > 1. 

For any / G L^( 


Since PfD(3{f) = 0 and P±D(i{f) = equation ( |7.61[ ) is equivalent to the 

system 

(7.62) i^i { - p + C\fi) - Sfi) m - s) - m-, d)} = 0, 

(7.63) Vm + P±{{&-V + C\yi) - d/r) ^(C)} = S)+M{^-, 6 


Our strategy will be to first solve (7.63) for /3 = d], for <5 > 0 and sufficiently 


small. We then substitute P[iJ,,S] into (7.62) to obtain a closed scalar equation. 
This equation is then solved for /i = /r(d) for 5 small. The first step in this strategy 
is accomplished in 
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Proposition 7.15. Fix M > 0. There exists i5o > 0 and a mapping G 

Rm.So = {|mI < X (Oi^o) which is Lipschitz in fi, such 


that l5{-\ pL,5) solves (7.63) for {pL,S) G Rm,So- Furthermore, we have the bound 






<6 \ 0 < J < (5o. 


The details of the proof of Proposition 7.15 are similar to those in proof of 
Proposition 6.10 in [^; equation (7.63) is expressed as (/ + p,,5) = 

V~^P± +A/’(^;(5)| and the operator C^{fx) is proved to be bounded on 

and of norm less than one for all 0 < 5 < i5o, with Jq sufficiently small. In 
bounding on require F{^{R) bounds for wave operators associated 

with the Dirac operator, V. These can derived from corresponding results for scalar 
Schrodinger operators, under the assumptions implied by n((') being a domain wall 
function in the sense of Definition 15.11 


7.6. Final reduction to an eqnation foru = and its solntion. Sub 


stituting the solution /3(^;/r,(5) (Proposition 7.15) into (7.62), yields the equation 
J'+[p,S] = 0, relating /i and S. Here, ff+[p,]S] is given by: 








L 2 (K) 
'Si 






The mapping {p,,5) G {l/rj < M, 6 G (0,5o)} is well defined and 

Lipschitz continuous with respect to fi. In the following proposition, we note that 
J+[p,S\ can be extended to a continuous function on the half-open interval [0,^o)- 

Proposition 7.16. Let do > 0 be as above. Define 


/or0<5<5o, 
[ fx- Po ford = 0, 


where po = — (a*, = E^‘^\ and is given in (6.27); see also (6.25) and 


(6.28). Fix M = max{2 |/io|, 1}. Then, {p, S) G {|/i| < M, 0 < S < <5o} t 77 {p, S) 
is well-defined and continuous. 


Proof: The proof parallels that of Proposition 6.16 of 12 . The key is to establish 


the following asymptotic relations, for all 0 < (5 < Jq with 6o sufficiently small: 


lim 6 (a*(-), Al(-; S) 
s^o \ ^ Vl2( 


= 1 : 


lim (5 ((J*(-),A/'(-;^) 

5^0 \ /L2(R) 


(7.64) 

(7.65) 


= -To', 
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and the following bounds hold for some constant Cm'- 


(7.66) 

(7.67) 

(7.68) 


-V) P{-,ix,5) 


(5(/i)/?(-;M><5) 


L 2 (K) 


L 2 (H) 


5 

\ / ( 


< CmS^-''-, 


< CmS''-, 

< Cm6. 


The detailed verification of (7.64)-(7.68) follows the approach taken in Appendix 
H of [l2]^_We make a few remarks on the calculations. Each of the expressions in 
(7.64)-(7.65) consists of inner products of the form: 


(7.69) 3(-5) = <5(a*(0,x(iei<5^"^) 


^{4>+(x,5^), J(x,(5i^2 •x))i2 
{4>_(x,(5C), J(x,(5i^2 •x ))^2 




Here, J(x,^) = e**^ ’^A(x, ^), where x i-A A(x, ^) is A^— periodic and ( i-A /C(x,^) 
is smooth and rapidly decaying on K. Consider, for example, the expression within: 


($+(x, •x ))^2 

_ ‘’ll 

Lemma 17.91 


This may be rewritten and expanded, using 


f P+{^,SO JC{^,Sfi 2 • x) dx = ^ / P+(x,dC) 1C (x, ^ +C) 

JS riGZ"'^ 

= f P+(x,<5C) ^(x,C) dx + ^ f p+(x,d^) A (^x, ^dx. 


Since in (7.69) ^ is localized to the set where |^| < d'" u > 0, for \n\ > 1, we 
have n/5 + ^ « n/5 and the decay of C i-t A(x, ^) can be used to show that, as 
5 tends to zero, the sum over \n\ > 1 tends to zero in L‘^{d^). It can also be 
shown, using the localization of that the n = 0 contribution to the sum, tends 
to ( P+(x), A (x, ^) ) . Therefore, uniformly in |^| < we have 


s^o 

Therefore, 


lim 6 ($+(x,d^), J(x,dl ^2 ■x ))^2 


= (^’+(x),/C(x,^) 


L 2 (n) 


limij;(d) = 

./b 


a*,+ (0 (-P+(x:),A(x,^) 


L 2 (n) 


(7.70) 


+a*-(C) (^’-(x), A(x,^) 


L2(n) 


The principle contribution to the limit in (7.64) comes from the A^i(^;d) term 
in (7.55). We apply ( |7.70 ) with the choice J = Jm{x-X) = and 

^Ai(x,C) = e"‘*^'’"JA4(x,C) = a*.+ (C)^+(x:) + a*-(C)P-(x) . 


The principle contribution to the limit in (7.65) comes from the Ai(^;d) and 
A/ 2 (C;d) terms in (7.56). We apply ( |7.70 l with the choice 

J = JAr(x,C) = (2il2 • Vx i9c - «(C)W^(x))i/'^^Hx,C) + 5 ^?/’^°Hx,C) 


dx 
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and C) = e C). The detailed computations are omitted since they 

are similar to those in [T^ . 

By E .64[ )-( |7.68| , it follows that f7+(^; (5) = fi — iiq + o(l) as 5 —)■ 0 uniformly for 
\fi\ < M. Therefore, J'[^J., 6 ] is well-defined on : \fj,\ < M, 0 < 6 < Jq}, 

continuous at 5 = 0 and Proposition |7.16| is proved. 


7.15 


to 


Summarizing, we have that given /3(-,p,, (5), constructed in Proposition^ _ 

complete our construction of a solution to ( |7.62 1-( 7.631, it suffices to solve ( |7.62 1 for 
fi = Furthermore, we have just shown that (7.62) holds if and only if ^ = fi{S) 
is a solution of i5] = 0. From Proposition 7.16 it follows that ^] = 0 has a 
transverse zero, ^ for all 6 and sufficiently small. The details are presented 

in Proposition 6.17 of [T^: 


Proposition 7.17. There exists Sq > 0, and a function 6 i—>■ /i((5), defined for 0 < 
(5 < (5o such that: |^((5)| < M, lim 5 _,,o /i(5) = /r(0) = ixq = and J[^{5),S\ = 0 
for all 0 < S < So ■ 


We have constructed a solution pair ^(5)^, with d^), of the 

band-limited Dirac system (7.52). Now apply Proposition 7.13 and the proof of 
Theorem |7.3| is complete. 


8. Edge states for weak potentials and the no-fold condition for 

THE ZIGZAG SLICE 


In Sectionj^we fixed an arbitrary edge, Oi = aiVi-|-a 2 V 2 and proved the existence 
of topologically protected Oi— edge states under the spectral no-fold condition. In 
this section, we consider the special case of the zigzag edge, corresponding to the 
choice: t)i = vi. We prove that the spectral no-fold condition holds in the weak 
potential regime, provided eVip > 0 ; this implies the existence of a topologically 
protected family zigzag edge states. 

We proceed in this section to prove the following: 


( 1 ) 

( 2 ) 

(3) 

(4) 


Theorem |8.2[ The operator — A -|- eV satisfies the no-fold condition along 
the zigza g (k 2 ) slice at the Dirac point ; see Definition 7.1 

Theorem |8.3l —A -|- £P(x) -|- dVF(x) acting in L^(Efc||=K.vi) has a spectral 
gap about the energy E = El. 

Theorem |8.4[ If eVi^i < 0, then the spectral no-fold condition for the zigzag 
slice does not hold. 

Theorem |8.5[ For 0 < |5| ^ and £ sufficiently small, the zigzag edge 
state eigenvalue problem for = —A -|- £E(x) -|- 6 K{S'k 2 ■ x)lF(x) has 

topologically protected edge states. 


We begin by stating our detailed assumptions on E(x) and bF(x). There exists 
Xo G such that E(x) = E(x — Xq) and lF(x) = VF(x — Xq) satisfy the following: 


( 8 . 1 ) 

(VI) 

(V2) 

(V3) 

(V4) 


Assumptions (V) 

Ah- periodicity: V(x -|- v) = V(x) for all v G A/j. 

Inversion symmetry: V(x) = V(—x). 

27r/3-rotational invariance: V(R*x) = V(x). 

Positivity of Fourier coefficient of eV, £Vip: £Vi,i > 0, 
where Vi_i = e“hki+k 2 ) y(7('y)^y. ggg (3 7 ) (2.3). 
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(8.2) Assumptions (W) 

(Wl) Ah- periodicity: + v) = iy(x) for all v € Ah- 

(W2) Anti-symmetry: W(—x) = —W{'s). 

(W3*) Uniform nondegeneracy of W : Let (x), j = 1,2 denote the ^, respec¬ 
tively, modes of the degenerate L^— eigenspace of = —A-\-eV. 

Then, there exists Oq > 0, independent of e, such that for all e sufficiently 
small: 


(8.3) 





> 00 > 0 . 


N.B. Consistent with our earlier convention, in the following discussion, we shall 
drop the “tildes” on both V and W. It will be understood that we have chosen 
coordinates with Xq = 0. 


Remark 8.1. We claim that (W3*) (see (8.3)) uniform non-degeneracy of W is 
equivalent to the assumption: 

(W3*) Wo,i + Wi,o - lUi,i ^ 0, 

where denote the Fourier coefficients of VF(x). To see this, note by 

Proposition 3.1 of [^, that for sufficiently small e, 

1 


ci>f(x) = 


iK 


VW\ 

see also (|3.5|). Evaluation of dj gives 




0 (e) 






lF(x) dx -I- 0(e) 


= 3 [(W'o.i + Wi.o) (r - r) + Wi,i(r2 - r^)] + 0(e) 

/ O 

= [W'o,! + lUi,o - + 0 {e), 


which is nonzero if (W3*) holds and e is sufficiently small. 

Let (K, E^) denote a Dirac point of = — A -|- £U(x), guaranteed to exist 

by Theorem |3.5| for all 0 < |e| < eo and assume that U(x) and VF(x) satisfy 
Assumptions (V), (W); see (8.1|-(8.2). 

In our next result, we verify the spectral no-fold condition for the zigzag edge. 
This is central to applying Theorem 7.3 to prove our result (Theorem |8.5[ ) on the 
existence of a family of zigzag edge states. 


Theorem 8.2. There exists a positive constant, 62 < Sq, such that for any 
0 < |e| < S 2 , = —A -I- eV satisfies the spectral no-fold condition at quasi¬ 

momentum K along the zigzag slice; see Definition \7.1\ 

By Assumptions (V), E_(A) < £’+(A) < Eb{X), b > 3. For any a sufficiently small: 


b = ±: 
b>3: 


a < |A| < 
|A| < 1/2 


1 


Ef{X)-El 


> 


|Li,i e] A > Cl |e| 


El’°{X)-El > C2 |£| (1 + | 6 |) . 
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Theorem |8.2| controls the zigzag slice of the band structure at K, globally and 
outside a neighborhood of Since a small perturbation of V which breaks 

inversion symmetry opens a gap, locally about K (see [13| ), it can be shown that 
H^s,5) g spectral gap about E = E^: 

Theorem 8.3. Assume T(x) and lT(x) satisfy Assumptions (V) and fWJ. Let 
be as in Theorem 8.2. Then, there exists C), > 0 such that for all 0 < |e| < £2 o,nd 
0 <6 <c\, e'^, the operator — A + eV {x) + SW (x.) has a non-trivial 
spectral gap about the energy E = E^. 


In the case where (V4) does not hold and the Fourier coefficient of eV, eVi,i, is 
negative: eVi^i < 0, then we prove that the no-fold condition does not hold: 

Theorem 8.4. [Conditions for non-existence of a zigzag spectral gap] Assume hy¬ 
potheses (V1)-(V3) but, instead of hypothesis (V4), assume eVip < 0. Then, for 
any e sufficiently small, the no-fold condition of Definition \7.1\ does not hold for 
the zigzag slice. 

The proofs of Theorems |8.2| and |8.3| are presented below. Section |8.1| discusses a 
reduction to the lowest three spectral bands. The general strategy, based on an 
analysis of a 3 x 3 determinant, is given in Section |8.2[ Theorem |8.4| is proved in 
Section |8.4| as a consequence of Proposition |8.15| 

We prove the following theorem on zigzag edge states using results of Section 

Theorem 8.5. Let = —A -|- eV{x.) -\- 5K{5\i.2 ■ x)IF(x), where P(x) and 

VF(x) satisfy Assumptions (V) and (W), and k(X) is a domain wall function in 
the sense of Definition \5.1\ Let £2 > 0 and c\, > 0 be as in Theorem \8.2\ and assume 
0 < jel < £2 and 0 < |i5| < ci,£^. Then, there exist edge states, 4>(x;fc||) € 
with |/c|| — 27r/3| sufficiently small. Furthermore, continuous superposition in 
yields wave-packets which are concentrated along the zigzag edge. 


Proof of Theorem \8.5\ We claim that the theorem is an immediate consequence of 
the spectral no-fold condition for —A -|- eV for the zigzag edge, stated in Theorem 


8.2 This follows by an applica tion of Theorem |7.3| (and Corollary 7.41. Since the 
details of the proof of Theorem 7.3 are carried out for the case of with e = 1, 

we wish to point out how the proof applies with e and 5 varying as in the statement 
of Theorem 18.51 


The proof of Theorem 7.3 uses a Lyapunov-Schmidt reduction strategy where 
the eigenvalue problem is reduced to an equivalent eigenvalue problem (nonlinear 
in the eigenvalue parameter, E) for the Floquet-Bloch spectral components of the 
bound states in a neighborhood of the Dirac point (K, E].). Stated for the relevant 
case of the zigzag edge, this reduction step requires the invertibility of an operator 
{I — Qs) acting on where Qs is defined in terms of Floquet-Bloch 

components in (7.29) for K 2 = k 2 . 

It suffices to show that the =k vi (^) norm of Qs is o(l) as J J, 0. From (7.31) 


in Remark |7.8[ the operator norm of Qs is bounded by e(5), given by: 

1^1 , l<^l < i 

uj{6’')ci{V) (1 + | 6 |)c 2 (P) -w((5-)ci(T) C 2 {V)' 

By Theorem 8.2 the no-fold condition holds with modulus a;(a) = and constants 


e(<5)< 


ci{y) = ciy|i4.ie| 


C2(T) = C2|£|. 
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< 




1 |5| 


Therefore, with a — S", 

ci6)< ^ + A . _ 

"" Ci{V) C2{V) q'^ci\Vis\ |e| C2 |e| 

By hypothesis, |i5| < c\,e^. Hence, 0 < c(i5) < + |e| —0 as e —0 if we choose 


u e (0,1/4). This completes the proof of Theorem 8.5 


□ 


To prove Theorems |8.2| - 18.4| we introduce a tool to localize the analysis about 
the lowest three bands. Throughout the analysis, without loss of generality, we take 
e > 0 and satisfy the cases eVip > 0 and eHi,i < 0 by varying the sign of Vip. 

8.1. Reduction to the lowest three bands. In this subsection we show, via a 
Lyapunov-Schmidt reduction argument, that the proofs of Theorems |8.2| and |8.3| 
can be reduced to the study of the lowest three spectral bands. To achieve this, we 
consider several parameter space regimes separately. 

We start by considering £ = (5 = A = 0. In this case, = (V -I- iK)^ 

has a triple eigenvalue, = |Kp, with corresponding 3— dimensional L^(]R^/A/,)— 
eigenspace spanned by the eigenfunctions Pa-, for a — 1 , t , r; see Section [XT] 


Next, we turn on e, keeping d = A = 0. From Theorem 3.5 there exists Cq > 0 
such that for e G = {—eo,eo) \ {0}, the operator (K) = — (V -I- -|- 

eV{x) has a double L^(M^/A/j) eigenvalue, E^. Let E^ = |Kp. The maps E^ 
and e ^ El are real analytic for e G with expansions ( |3.10 ), (3.11): 

El =El + e(yop-Vi,^) + 0{e^), 

El = EI -|- e(Vo,o + 2Hip) -I- 0{e^). 

More generally, we may study the eigenvalue problem 
(8.4) - E)p = 0, pGl^iR^Ah), with E = El + p, 


and seek, via Lyapunov-Schmidt reduction, to localize (8.4) about the three lowest 
zigzag slices. Written out, the eigenvalue problem has the form: 

(8.5) \ - {V + i[K + Xk2]f - El - p + eV + SW ]p{^) = 0, p G L'^ / A^). 

Since e and 6 will be chosen to be small, we shall expand p(x) relative to 
the natural basis of the L^(M^/Aft,)— eigenspace of the free operator, = 

— (V -I- iK)^, displayed explicitly in (3.5). Let pH denote the projection onto 
span {pa- : a = and P-*- = I — PK 

We seek a solution of (8.5) in the form p(x) = p||(x) -|-p_l(x), where 

r 1 

P|| e spanjpcr : cr = l,r,r} , p± G span : o’= 1, t, r} 

Then, we have that (|8.5|) is equivalent to the coupled system of equations: 


( 8 . 6 ) 


-{V + i[K + Xk2]y -El-p 


P\\ 


+ 


ePW + SP^^W 


P\\ 


eP^^V + SP^^W 


P± = 0, 


- {W + i[K + Xk2]f - El- p 


P± 




eP^V + 5 P^W 


P± 


eP-^V + SP^W 


= 0 . 


(8.7) 
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Proposition 8.6. There exists a constant c > 0 such that if jej -|- |^| < c then 


= r_(v+i[K + Ak2])"-£;:-Ai 


is invertible on P'^L^l 


VA.). 


Proof of Proposition\8.()[ This follows from the lower bound: 


0 < On = inf inf 

|A|<l/2 m^{(0,0).(0.1),(-1.0)} 

where mk = ruiki + m 2 k 2 . 


|K + mk + Ak2p- IKl 


□ 


By Proposition 8.6, for all e and /x sufficiently small equation (8.7) is equivalent 


to: 


I + 


^(e.o.A) _ 1 [spJ-y + 5P-^W] 


P± 


p(e.o.A) _ J [eP^v + 6 P^W] p||. 


Suppose now |e| +15| + |/x| < di, where 0 < di < c is chosen sufficiently small. Then 
we have 


P± = 


( 8 . 8 ) 


p-(£,o.A) _ J [ep-Ly + SP-^W] 


p-(^,o.A) _ 1 [ep-Ly + SP^W] P|| = P{e, S, A, fi) 


P\\- 


Equation (8.8) defines a bounded linear mapping on P^^L'^(Mf/Ah) with operator 
norm < 1: 

P\\ ^ P_L[e, A;p||] = V{e,6,X,fj.) py : Ran Ran (P-^) , 

which is analytic in e, S, A and p for |e| + |(5| + |/x| < di and |A| < 1/2. Consequently, 
equation (8.6) becomes a closed equation for py which we write as: 

M{e,S, X, pL)p = 0, 

where 

M{e, S, A, ^l) = Pll f - (V + i[K -h Ak 2 ])^ -El-fi + sV + 6W 

pll 


(8.9) 


( 8 . 10 ) 


+ {eV + SW)P{£, (5, A, fa) 


= pll 


- (V + i[K -f Ak2])^ -El-fi + eV + SW 


+ {eV + SW)P-^ [pC-oP) - f^ + sV + SW] P-^{eV + SW) 


pll 


The operator Ai{e,S, X, fi) acts on the three-dimensional space pllp^(M^/A;i) = 
spanjpo- : a = l,r, r}. Moreover, from the expression (8.10) it is clear that for 
e,6,X,fi complex, M{£,6, X, fi) is self-adjoint. We shall study it via its matrix 
representation relative to the basis {pa- : cr = 1,t, r}: 

(8.11) M^^a{£, 6 ,X,fj.) = {pa,M{£, 6 ,X,fa)pa)i^ 2 (j^y^^), a = l,T,T. 


Clearly, M(e, S, A, p) is Hermitian for e, <5, A, fi real: 5, A, fa) = S, A, fa). 

Summarizing the above discussion we have the following: 

Proposition 8.7. (1) The entries of the 3x3 matrix M (e, 6, A, fo) are analytic 

functions of complex e, 5, A and fa for |e| -|- |(5| -|- \fa\ < di and |A| < 1/2. 
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(2) For jel + |(5| + \fi\ < di and |A| < 1/2, we have that E = is an 

/Ah)— eigenvalue of if and only if det M{e,d, X, fi) =0. 


We now study the roots of det M{e,S, X, fj,) = 0 (eigenvalues of for e 

ae 

), 


and S small. First let e = <5 = 0. By the formulae ( |8.42 1 and ( |8.46 ), derived and 
also applied in Section [STSl we have: 


M{0,0,X,fi) = (^ {p^ , M{0,0,X,fj.)pa) 


L 2 (R 2 /A;,) 


a.a—l,T,T 


( (Pa, (-(V + *(K + Ak2)f),p,,)^,(^,/^^) - (F;/ + M)(5..a) 


(T^a—l 


' — p aX aX 

aX X^q^ — p aX 
aX d-A X^q^ ~ Pj 


a = IT. 

73 


Thus, M(0,0, A, p) is singular if /x is a root of the polynomial: 


det M(0,0, A, p) 

= -7 + + pidX^W]^ - 3A^7) + - 3A'‘gVl^ + 2A33?(a7 

= -p^ + 7(3A^7) + - 3A"‘7) + A®7 + 


where we have used that \a\^ = ^ and Sfi(a^) = 0. The roots, p^^\x), j = 1,2,3, 

defined by the ordering p^i\x) < p^\X) < p^^\x), listed with multiplicity, are 
given by: 

(8.12) 0 < A < i : p^°\x) = 7A(A - 1), Ai2°^(A) = g^A^ pf\x) = g^A(A + 1), 

(8.13) -^ < X < 0 : p^°\x) = g^A(A + 1), ^2°^^) = g^A^ pf\x) = g^A(A - 1). 

The roots p^^\x),j = 1,2,3, are eigenvalues of — (V + f [K + Ak 2 ])^ — i?/. They 
are uniformly bounded away from all other /Ah)— spectrum. This distance 

is given by 


go = inf 

|A|<l/2 


inf 

ri^{(0,0),(0.1),(-l,0)} 


|K + mk + Ak2r- |K| 


>0, (i?° = |K|7, 


where mk = ruiki + m 2 k 2 . 

Note that det M{0,0,0, p) = —p^ has a root of multiplicity three: p = 0. By 
Proposition 8.7 and analyticity of det M(e, 6 , A, p), we have that for e, S and A in a 
small neighborhood of the origin in C^, there are three solutions of det M (e, 5, A, p) = 
0. We label them pj{e,6,X), j = 1,2,3. By self-adjointness of M{s,S, X, p), these 
roots are real for real values of e and S. Correspondingly, for small and real e, <5 
and A there are three real ^fc||= 27 r/ 3 ~ eigenvalues of denoted 

eI^'^\x) = Ef’^\K + Ak 2 ) = El + p,{e,6, A), j = 1,2,3. 

The ordering of the Ej implies 

Piie,6,X) < p 2 ie,S,X) < p3{e,S,X). 

A mild extension of Proposition |4.5| yields 
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Proposition 8.8. For each |A| < 1/2, there exist orthonormal eigenpairs 

($± (x; A), i?±(A)) and (<i>(x; A), i?(A)), real analytic in A, such that 

span {<i)_(x; A), ‘&-(-(x; A), ‘&(x; A)} = span {<i)_;(x; K -|- Ak 2 ) : j = 1, 2,3} . 

For fixed e small and A tending to 0 we have: 

( 1 ) 

(8.14) Ef=\X) = El ± |Af| |Z 2 | A + A^ e±(A,e), 

(8.15) = El + A2e(A,e), 

where e±(A,e), e{X,e) = 0(1) as X,e ^ 0, and Z 2 = + ik^\ These 

functions can be represented in a convergent power series in e and X in 
a fixed neighborhood of the origin. Furthermore, e±(A,e), e{X,e) are 
real-valued for real X and e. 

(2) Therefore, for all small e and X, the three roots o/det M(e, <5 = 0, A, /j,) may 
he labeled: 


(8.16) /z+(A,e) = |A|| \z 2 \ A + A^ e+(A,e) 

(8.17) /i_(A,e) = —|A|| \z 2 \ A + A^ e_(A,e) 

(8.18) fi{X,e) = El — El A^ e(A,e) , where — i?* = SeVi^ + O(e^). 


Proof of Proposition \8.S^ Part 1 can be proved as follows. The expansion (8.14) and 
analyticity follow by perturbation theory as in Proposition 4.5 see also |15[|22| . The 


expansion (8.18) also follows by perturbation theory of the simple eigenvalue El 
for A = 0. It is easy to see that 

(8.19) i?'^’^=°(A) = El + X X -2*k2-^ + A^e(£,A). 

We claim that (Vx4>®\ = 0. This follows since <i>® G ^ as in the proof of 

\ / Z/T^ ’ 


Proposition 4.1 of 13 . This proves the expansion (8.15). 


□ 


Finally we discuss a useful symmetry of det M(e, 6 , X, p. = 0). 


Proposition 8.9. Assume V satisfies Assumptions (V) and W satisfies Assump¬ 
tions (W). Recall M{e,S,X,0), defined in (8.9)-(8.11). Then, det M(£, 5, A, 0) is 
real-valued for real e, 6, X and analytic in a small neighborhood of the origin in C^. 
Furthermore, 


(8.20) det M{e, 5, X, 0) = det M(e, -5, A, 0) 

and therefore det M(£, 5, A, 0) is a real analytic in e, S'^ and X, and we write 

D{s, S‘^,X) = det M(e, d, X, 0) 


Proof of Proposition \8.9[ Let I[/](x) = /(—x) and C[/](x) = /(x). Using that e, 5 
and A are real, V (x) is even and W (x) is odd, one can check directly that 

(8.21) Colo oXoC 

Furthermore, note that 

(8.22) (Co I)pi =pi, (C oI)p^ = p^, (C oI)p^ = p^. 
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It follows that 
(8.23) 


Colopll = pll oCoJ. 


Using the symmetry relations (8.21)-(8.23) to rewrite M(e, —<5, A,0) in terms of 
we find that M{e,—6,X,0) can be transformed into M{e,d, X,0) by inter¬ 
changing its second and third rows, and then interchanging its second and third 
columns. Therefore, det M(e, d. A, 0) = det M(£, —<5, A, 0) and the proof of Proposi¬ 
tion is complete. □ 

8.2. Strategy for analysis of det M(e, d. A,/r) in the case eVi^ > 0. We first 
observe that for a positive constant, di, if |e| -I- |5| < di then 

(8.24) \E^^’^\X) - El\ >C4(l + |j|), j>i. 

Indeed, by the discussion following Proposition |8.7| we have that there exists d 2 > 0 
such that for j > 4, A)| = |uj^’‘^^(A) — > ^ 2 ; the lower bound (8.24) now 

follows from the Weyl asymptotics for eigenvalues of second order elliptic operators 
in two space dimensions. 

Hence, we restrict our attention to Ej‘^’^\x) = E^ + j = 1,2,3, which 

we study by a detailed analysis of det M(e, i5. A,/i = 0). The analysis consists of 
verifying two steps, which we now outline. 

Step 1: Fix C|, >0 and arbitrary. We will prove that there exists > 0, such that 
the following holds. There exists ei > 0 and constant C 3 , depending on V and W, 
such that for all 0 < |e| < ei and 0 < |5| < c\, e^: 


(8.25) ^ 1^1 ^ 2 


\E, 


{e,S) 


iX)-El\ >C3e, j = 1,2,3. 


Furthermore, by ( 8.25| ) and (8.24), it follows that 
(8.26) 


C\;^/e < |A| < - 


L^(M/A,,)-spec(i7('^’^’^)) n 


El — ce, El -I- ce 


is empty. 


Step 2: Let and C\, be as in Step 1. We will prove that there exists 0 < £2 < £ 1 , 
such that if (A, S) are in the set: 

(8.27) |A| < Cb £5 and 0 < |i5| < c\, where 0 < |£| < £ 2 , then 

(8.28) 

det M{s, A, (5,0) = (^^A^ -I- £Vi^i^ 6 ^ |lFop -I- VFi^o ~ Wh,i| ^ ^ ( 1 + o(l) ). 

A simple lower bound on the three eigenvalues |/rj(£, A)| = |A'j(A) — El\ is then 
obtained as follows. By (8.28), for some positive constants: Ci and C 2 , we have 

Ci(A^-te) • (A^< |detM(£,A,(5,0)| 

= |det M(e, A, 5, 0) — det M(£, A, <5, Pj (A, £, (i))| < C 2 \0 — Pj (£, A A)| = C 2 |Pi(£, 5, A)|, 
for j = 1, 2, 3. Therefore, with C 3 = CxjC^-, 

(8.29) Ef{X)-El =\^x,{e,5,X)\>C^e{X^ +5"^) > C 3 e j = 1,2,3. 


It follows from (8.29) and (8.24) that 

the L‘^(M.‘^/Ah)— spec( ^ jg empty 

with rj = 4(73 £ whenever £, {X,S) satisfy the constraints (8.27). 
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Theorem 8.2 and Theorem 8.3 are immediate consequences of (8.26) (Step 1) 
and (8.29) (Step 2) and the representation: 






rr/a) 




'\M<h 


L2(BVAh) 


dX. 


Hence, we now turn to their proofs. We first carry out Step 1 by a simple 
perturbation analysis about the free Hamiltonian, We then turn to Step 

2, which is much more involved. 

Verification of (8.25) from Step 1. Let C\, denote a positive constant, which 
we will specify shortly, and consider the range C\,^ < A < 1/2. Using the expres¬ 
sions for ^j°^(A), j = 1,2,3, in (8.12) we have that if £ < e'{Ci,) = then 

\e[°Hx) - E°\ = q^\X\ |1 - A| > C,qWe/2, juf (A) - U/| = q^\X\^ > C^q^e, 
and \E^\x) — E^\ = g^|A| |1 -I- A| > C^q‘^^/ej2. Note that the eigenvalues 
(£,(5, A) i-A- E^’^{X) are Lipschitz continuous functions; see Chapter XH of or 
Appendix A of [Id]. Therefore, we have 


\E. 


t^\x) - El\ >C^q‘^e- |£| 


- 1^1 \\W\\^ > - q^ e, 


for some C\, positive, finite and sufficiently large. With this choice of C\,, we also 
have 

\e[^’^\x) - El\ >a V^/2 - |e| 

.CA)(A)_^e| - |£| 


\Ek 


- H IllViloo, 

- H ||HA|U. 


Therefore, there exists £i > 0, such that for all 0 < |£| < £i we have 
\e[^’^\x) - El\ + \Ei^’^\x) - El\ > C,q^^e/4. 

This completes the proof of the assertions in Step 1. 

8.3. Expansion of M{e,5,X,0) and its determinant for £Vi.i ^ 0. The key to 
verifying Step 2 and proving Theorem |8.2|is the following: 


Proposition 8.10. Let C\, be as chosen in Step 1; see (8.26). Then, there exist 
constants £2 > 0 and c > 0 such that for all 0 < e < 62 , if 

(8.30) 


(8.31) 
where 

(8.32) 


0 < |A| < C\, £^/^ and 0 < |(5| < c £^, then 
— det M(£, (5, A, 0) =7r(£, i5^,A) -I- o ((A^-|-£)(A^-I-<5^)) , 

^{e,6\X) = (g2A2+£Vi,i) (q^X^ + 5^ jw'o,! + Wyo - lUi,i|") ■ 


Elere, Wm, m G if', denote Fourier coefficients ofW{x) and, by (8.30), the cor¬ 
rection term in (8.31) divided by (A^ -I- £)(A^ -I- i5^) tends to zero as e tends to zero. 
Thus, 

£Ui 4 > 0 — det M{e,S, X,0) =tt{£,S‘^,X) (l-l-o(l)) in the region { 8 .SO). 

We now embark on an expansion of det M(£, (5, A, 0) and the proof of Proposition 


8.10 M{e, S, A, 0), see ( |8.9[ )-(8.11 ), may be written as the sum of matrices: 
M(£, S, a, 0) = [M° + M^ + ] (£, S, A, 0), 


(8.33) 
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where the a,a = 1, r, r entries are given by: 

(8.34) MO^(e,5,A,0) = /p,,(-{V + t[K + Xk2]f-El)p^ 

\ V / /L2(K2/Afc) 

^ (PcTT^P^) L^(R^/Ah) ^ 

^ {Pcri^P^) L^(R^/A,^) 1 

= (P<T,(e‘^^ + ^W^) P(e,(5,A,0) ps)^2(r2/a^) . 

For e and 6 small, 

(8.35) M°^(e,^,A,0) = M°;f^™"(e,A) + 0(s^), 

where (inner products over L^(R^/A/i)) 

(8.36) 

M“;7^”"°"(£, A) = (p., (- (V + i[K + Aka])" - [F;*“ + e{Vo,o - Fi.i)]) ps) , and 

(8.37) 

M^j(£,5,A,0) 

= (^{eV + 5W)p^,P^ {-{V+ i[K + \k 2 ]f - El + eV+ SWy^ P^{eV+ SW)pa') 
= (^{eV + 5W)pa,P^ {-{V + i[K + \k2]f - Ety^ P^ieV + 5W)pa) 

+ 0(5^ + 5^e + fc" + e®) = 0(e" + ed + 5"). 


We next explain that to calculate the determinant of M{e,6,X,0) to the desired 
order in the region (8.27), it suffices to calculate the determinant of the approximate 
matrix: 


(8.38) M“PP’'“(e,5,A) = M°’“pp™^(£, A) + M'^(e) + M’^((5). 

That is, we show that the omitted terms in M{e,S, X,0) — M‘^pp^°^{s,S, X,0) con¬ 
tribute negligibly to the determinant of M(e, i5. A, 0), when compared with the poly¬ 
nomial, 7r(e,5^, A), in (8.32), provided A and S are in the region ( |8.27| ): 

(8.39) |A| < C[, £2 and |(5| < ct e". 


where C\, and ci, are appropriately chosen constants. 

Recall that D{e,6^ = 0, A = 0) = det M(£, 0,0,0) = 0 for all e, since p = 0 
corresponds to £1 = El, which is an eigenvalue of 7 j(e-^=OA=o) _ Thus, 

D{e,5‘^,X) = det M(e, (5, A, 0) is a convergent power series in £, 5'^ and A with no 
“pure £” terms. On the other hand, the entries of the matrix M{e,6, X,0) are 
convergent power series in e, S and A. 


Proposition 8.11. For (A, i5) in the region (8.39) we have 

|A| 6 ^ = o{{X‘^+e)iX‘^ + S^)). 

Therefore we may drop the 0{XS'^) terms, a further simplification. 


Proof of Proposition \8Jl\ Consider separately the two regimes: (a) |A| < and 
(b) |A| > £^-^. For |A| < e^'^, we have |A|5^ < e^-^5'^ = £°-^eA^ < £°-^(A" -|-e) • (A" -I- 
(5"). And for |A| > e^’^, note that (A" -I-e) • (A" -I-(5^) > eA" > £• e^-^ = e^-". On the 
other hand, |A|J^ < < £^ < £’® • (A^ -I- e) • (A^ -I- <5^). This completes the proof 

of Proposition [8TT] □ 
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Let US now attach “weight” 1 to the variable e and “weight” 1 /2 to the variables <5 
and A. A monomial of the form carries the weight a+&/2-|-c/2, where a,b,c G 


N. In Proposition 8.12 we show that all terms in the power series of det M(£, A, (5) 


which are of weight strictly larger than two introduce negligible corrections to 


det A, (5) for A and J in the region (8.27|. 


Recalling that there are no pure £ terms, we see that a monomial in the power 
series of D{e, S^, A) of weight larger than 2 must have one of the following two forms 
(a, b,c G N): 

A X 

X e‘^X\S^y, 


(I) 

(II) 


with a + 6/2-|-c > 3/2 
with a + bl2 -|- c > 1 


2a + 6 + 2c > 4 
2a A 6 A 2c A 3. 


Proposition 8.12. Terms of form (I) and form (II) that may appear in D{e, <5^, A) 
are o ((A^ A £)(A^ A 6"^)) as £ — )■ 0, for (A, i5) in the region (8.39); |A| < C\, £5 and |<5| < 


£^. We may therefore neglect all terms in the power series of D(e, , X) which are 


of weight strietly larger than 2 for (A, in the region (8.39). 


We prove Proposition 8.12 by estimating all terms of the form (I) or (II), and we 
begin with the following lemma, which is a consequence of part 2 of Proposition 


Lemma 8.13. Let/i_(£,A), fj,+ {e,X) andJl{e,X) denote the three roots of det M {e , S = 
0, A,/i), defined and analytic in (£, A) in a neighborhood of the origin and dis¬ 
played in (8.16)-(8.18). Then, 


det M(£, 6 = 0, X, fj.) = {fi — X)) X {p — p,+ {e, X)) x (p — A)) x Q{e, A, p.), 

where n{e,X,p) is bounded. In particular, for all e such that 0 < |£| < £o there 
exists a constant Cr such that 


(8.40) 


|det M(£, (5 = 0, a, 0)1 < Ce |A|^ 


Proof: For hxed £ and A, the mapping p i-A det M(£, (5 = 0, A, p) is analytic for |^| < 
po with zeros at p±{e, A), p{e, A); see Proposition |8.8| . Fix e' and A' small such that 
if |£| < e' and |A| < A', the roots all satisfy |^+(£, A)|, |/i_(£, A)|, |/I(£, A)| < po/^. 
We claim that for such e and S, 

_ detM{e,S = 0,X,p) _ 

~ (M-Ai-(e,A)) X (p-p+{e,X)) X {p-p{e,X)) 

is uniformly bounded for all |^| < po, |£| < s' and |A| < A'. Indeed, since 
the roots are bounded in magnitude by po/2, we have max|^|^^Q |n(£, A,/i)| < 
(2/pQy maxj^i^^g \det M{e,6 = 0, A,/i)|. Applying the maximum principle we have 

max |r2(£, A,^)| < (2/^o)^ max |det M(£, 6 = 0, A,/r)| < {2/poyC{po,e',X'), 
AI</jo \r\=uo 


where C{po, s', A') is a constant. The bound (8.40) now follows from the expansions 
of the roots. 

Proof of Proposition |8.12[ 

(I) Terms of the form A x £“A^(6^)“, with 2a A 6 A 2c > 4, a, 6, c € N: 

(i) Sup pose f irst that c = 0. Then, we consider A x £“A^ with 2a A 6 > 4. By 
Lemma 8.13 we must have 6 > 1. Thus, A x £“A^ = X'^e°'X^~^. If a > 2, then 
A X £“A^ = A^£^£““^A^“^ < (A^ A S"^) X (A^ A £)^ = o ((A^ A e)(A^ A <5^)) for (A, <5) 
in the region (8.39). Otherwise, a = 0 or a = 1. If a = 0, then 6 > 4 and 
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A X = A • A^ X A^ • A''-^ < A • (A^ + e) x (A^ + S‘^) = o ((A^ + e)(A2 + S‘^)) 

for (A, (5) in the region (8.27). Now suppose a = 1. Then, b > 2 and we have 
A X e“A^ = A • A^ X £ • A°“^ < A • (A^ + <5^) x (A^ + e) = o ((A^ + £)(A^ + 6"^)) for 
(A,d) in the region ( 8.39[ ). 

(ii) Suppose now that c > 1. If 6 = 0, then A x £“A^(d^)'^ = A(5^£“((5^)'^“^ with 


2a + 2c > 4, which is = o ((A^ + e)(A^ - 
8.11 Finally, if 6 > 1, then Axe“A*'((5^) 
for (A, S) in the region ( |8.39[ ). 


S^)) for (A, d) in the region by Proposition 
= A(52.e“A*'(A2)^-i = o ((A^ + £)(A^ + 6^)) 


(II) Terms of the form x £“A^(d^)‘^, with 2a + 6 + 2c > 3, a, 6, c S N: 

(i) Suppose first that a = 0. Then, S'^ x £°‘X^{S'^Y = S'^ x A*’((5^)°, 6 + 2c > 3. If 
& > 1, then we rewrite this as d^A x A^“^((5^)‘^ = o ((A^ + £)(A^ + d^)) for (A, <5) in 
the region ( |8.39 l, by Proposition 8.11 And if & = 0, then x £“A^(d^)'^ = (<5^)°+^, 


with c > 2, which is < < S'^£-£^{S'^y~^ < (A^ + d^)(A^ +£) •£^(<5^)°“^ = 

o ((A^ + £)(A^ + (5^)) for (A, 5) in the region ( |8.39| ). 

(ii) If now a > 1, then <5^ x £“A^((5^)^ = 6‘^e ■ £““^A^(d^)'^ < (A^ + (5^) • (A^ + £) • 
£a-iAf)(j 2 )c _ o ^(A^ _|_£)(a 2 + d^)) for (A,d) in the region (8.39). This completes 
the proof of Proposition |8.12| 

Now each entry in the 3x3 matrix, M(£,A, A,0) is a sum of terms of weight 
> 1/2; this is a consequence of the expansion (8.33), (8.34), ( |8.35 1, ( |8.8| ) and the 
explicit expansion of displayed in (8.47). If we change any one entry by a 

term of weight > 3/2, then the effect on the 3x3 determinant D{£, S^, A) will be a 
sum of terms of weight > 3/2 + 1/2 + 1/2 > 2. By Propositions 8.11 and 8.12 such 
terms are o ((A^ + £)(A^ + d^)) in the region (8.39). Therefore, we may compute 
each entry of M{£, 6, A, 0), retaining only terms of weight strictly smaller than 3/2 
and discarding the rest. The resulting determinant will differ from T)(£, i5^,A) by 
terms which are o ((A^ + £)(A^ + d^)) in the region ( 8.39[ ). 

We next study the power series of the 3x3 matrix, M (£, d. A, 0), keeping in mind 
that the relevant monomials are those of weight >1/2 but strictly less than 3/2. 
The complete list of such monomials is: £, d. A, A^, 5^ and Ad. 

Before proceeding further we show that in fact that a monomial of type 5^ 
can be neglected. Indeed, the weight < 2 contributions of such a monomial to 
D{£, a, 0) = det M(£, d. A, 0) will be a sum of monomials of the type: (i) d^ x d-d, 
(ii) d^ X A • A and (hi) d^ x A • d. Terms of type (ii) and (hi) are clearly o( A d^) = 


8.11 


o ((A^ + £)(A^ + d^)) as £ —>■ 0 for (A, d) in the region (8.39), by Proposition 
For the type (i) term we have, for (A, d) in the region (8.39), d^ x d • d < d^ £ £d < 
(A^ + d^) • (A^ + £) £d = o ((A^ + £)(A^ + d^)) as £ —>■ 0. Hence, we may strike d^ 
from our list. In particular, we may neglect the contribution from the matrix A4^; 
see (8.37). 

Relevant monomials in the expansion of M(£,d, A, 0) : We shall call the 
monomials: £, d. A, A^ and Ad relevant . All others are called irrelevant . 

Stepping back, we have shown above that M = + M'^ ((8.33l), 

where + 0{£^) (( |83^ ) and = 0{£^ + £d + d^) ( |o7| ). 

We have further shown that the relevant monomial contributions for calculation of 
det M(£, d. A, 0) are all contained in M“^'P’’°“(£, d. A, 0) = M°’“^’^"’°“(£, A) + M^(£) + 
M'^(d). We consider each of these matrices individually, and explicitly extract the 
relevant terms in each; see Proposition 8.14| below. 
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Expansion of ^o,approx, entries of are 


M 


0,approx 


{e,S,X,0) 


{pa, ( - (V + z[K + \k 2 ]f - [E° + e{Vo,o - ) Pa) 


A"(KVAfe) 


(8.41) =(p^,-{V + i[K + Xk2]fpa) - [E° + eiVop - Vi^i)] , 

\ /L2(R2/Ah) 


where we have used that {pcr,Pa) = S< 7 ,a- The first term in (8.41) may be written, 
using (3.5)-(|3.6|), E° = |Kp and |k 2 p = as: 


(8.42) {p^,-{V + i[K + Xk2]fpa) 

\ /L2(RVAfe) 

= (p<T,-(V +*K)^pa) - 2iAk2 • (po, (V+ *K)pa) + {pa,X‘^q^pa) 
~ ^<y,a + J(j,a- 


Consider now the matrix 


Ja,a = - 2 ik 2 • ($< 7 , V4>ar)^^ = - 2 ik 2 • J $o- V4>adx 

(8.43) = 2 k2 • i (/ + aaR + aaR^) K. 

We pause to collect some properties that will enable the evaluation of Ja,a', see 
also j^. Recall that R has eigenpairs: (t, ^) and (r, ^), where ( = ^(Ij*)^- 

Then, tR has eigenpairs: [1, C] and [r, ^]. Furthermore, tR has eigenpairs [1, and 
[t,C], and 

- [ / + ri? -I- (ri?)^] C = C) [ -^ + C = 0) 

i [/ + TR+(rR)2]C = C, [/ + ri?+(Ti?)2]C = 0. 

Hence, | [ f + tR + (tR)^] and | [ f + tR + (ri?)^] are, respectively, projections 
onto span{ ( } and span{ ( }. For any w G C^, we have w = {(, w)j ,2 C + (Cj Cj 
where (x, y)j ,2 = x • y. Therefore 

(8.44) i [/-hrR-f (ri?)2] w = (C,w)C, ^ [ / + rR + (ri?)^] w = (C,w) C 

Also, {I — R){I + R + R^) = I — R^ = 0 and therefore 


(8.45) 


I + R + R^ = 0. 


We next calculate Ja,a using (8.44)-(8.45). Note that J is Hermitian, and by (8.45) 
its diagonal elements Jo-,a all vanish: J^-^a = Ja,a, Ja,a = 0, cr = 1, r, t . It suffices 
therefore to compute the three entries Ji,r, Ji.r and Jry- 


Ji,r = 2 i [/ + TR+(Ti?)"]K.k 2 =2 (C-K) (C-k 2 )= a, 

Ti.r = 2 i [/ + TR+(ri?)"]K.k2 =2 (C-K) (^ k2) = 5, 

Jr,r = 2 i [/ + TR+(Ti?)2]K.k2 =2 (C-K) (C • k2) = O. 
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Thus, 


/O a a\ 

(8.46) J = ( Jcr,a- ) = I S 0 a j , where a = 2 (^ • K) ((^ • k 2 ). 

ya a 0 J 

It follows that 

(8.47) <5, A, 0) = ( -e(foo.o - ^i.i) + ) <5, ~ + A J^^a- 

Expansion of {s): M^-{£) =e {pa,V Pa)i2(R2/A^) = e V^^a, where 

Vo-.a- = (Pcr,Vpa}j^2(^g2/j\^-) 

= ^ {1 + era+ aa)Vop + aVop + aVo-i + aVi^ 


+ crfo-in + aaVii + aaV-i^i 


Since V is real-valued and even, it follows that V-m = Vm- Furthermore, V is also 
R— invariant and therefore Vb.i = ^i,o = ^i,i- Hence 

V(T,a =-(1-I-cr ft-I-CT cr) Vb,o + -(u-I-u-I-d-I-(f-I-(T(t-I-CT ff) Vi^. 
o o 

V is clearly symmetric and using that = l-|-r-|-T = 0, we obtain 

[e] = e V, where 

/Ho,o + 2 Vi^i 0 0 

V = 0 0 

V 0 0 Vo,o-Hiu, 

Expansion of (5): M^(e) = 6 {pa,W Pa)L^^m/A^) = ^ where 

(8.48) yVa,a = ^ aWoA + dWo,_i -I- ctIF_i^o+ o'IFi^o + crcHFi^ -I- . 

O L 

Since W is real and odd, we have that = —fFm and Wm is purely imaginary. 
Therefore, 

Wcr,a = A (o' “ ^0,1 + (o’ - ct) VFyo -I- (cr O’ - cr d) IFyi 


, a, a = l,r, r. 


It follows that M^{d) = 6 W, where 


0 


r —r' 


0 T —T\ 


0 -1 1 


w = Wol T —1 0 -f WlO T —1 0 -t-Wll —1 1 


0 


0 1 


0 1 


0 - 1 ; 


and Wij = —i Wij/y/i G M. 

Now assembling all relevant terms (weights >1/2 and less than 3/2) we obtain 

Proposition 8.14. For a,a = 1,t, r, 

M,^a{e, <5, A, 0) « ™"(£, <5, A, 0) 

« ( ~£(Fo,o ~ Hi^i) + }?) Sc^a + A J^^a + £ ^tT,a + ^ FVo-,a- 

Here, Aa-^a ~ Ho-,a, means that their difference is a matrix with entries having 
weight >3/2. Hence, the contribution of such terms to the determinant consists of 
terms of weight strictly larger than 2, for (A, (5) in the region (8.39). Hence, these 
terms can be neglected, by Proposition\8. 1^ 
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So the calculation of det M{s,S, X,0) boils down to the calculation of det (e, (5, A, 0). 

Calculation of det 5, A, 0): Assembling the above computations, we 

have that 

/A^g^-I-Seld,! a\ — 5w aX + 5w \ 

J^approx _ I aX—Sw A^g^ — (5(ui01 + Wio — Wll) qA j , 

\ aA-I-aA A^g^-I-(5(woi + wio — wii)/ 


where w = wn — wqit — wiqt and a = 2(C • K) (C • k 2 ). Note that: 

2 2 

(8.49) a = IT, 3fi(a) = —%r, 5i(Qf^) = 0. 

V 3 2 


Calculating the determinant of , and using (8.491 and that Wij = —iWij /-s/S 

yields: 


det <5, A,0) = - (g^A^ + eVi,i) (g^A^ + 3(5"(zi;2^ + wlo + «^ii)) 

+ 6eVi,ii5^(wiiWoi + wiQWii - woiwio) + 0{X5‘^) + 0{eX‘^) + 0{X^) 
= — {q^X? + eVip) ( g^A^ + 35^(woi + wio ~ wii)'^ ) 

+ 0{X^S^) + 0{XS‘^) + 0{eX'^) + 0{X^) , 

= - (g^A^ + (g^A^ + jR/o,! + Vbi,o - VCi,i |') 

+ 0{X‘^d‘^) + 0iXS‘^) + 0ieX^) + 0{X^) 

(8.50) = -7r(e,(5^,A) + o ((A^ + e)(A^ + 5^)) , 

for (A, S) in the region (8.39). This completes the proof of Proposition 8.10 


8.4. If eCi,i < 0, the zigzag slice does not satisfy the no-fold condition. 

Recall that to satisfy the case eVip < 0 we assume, without loss of generality, that 
e > 0 and Vi.i < 0. Theorem 8.4 follows from: 

Proposition 8.15. Assume 

(8.51) 0 < |£| < £ 2 , and 0 < <5 < ci, £^. 

There exists Oq > 0 and Ag > 0 satisfying e < X^ < 0o^/e such that for all e 
sufficiently small, 

det M(£, (5, Ae, ^ = 0 ) = 0 . 

Thus, Ef_ is an interior point of the (S)— spectrum of 


It follows that for £Vi,i < 0, the operator does not have a spectral gap about 

E = El along the zigzag slice. Referring to the middle panel of Figure we see 
that, for (5^0 and small, a local in X gap opens, for A small, about the energy 
E = El. But since the no-fold property is not satisfied (by Proposition |8.15[ ), this 
is not a true (global in A S [—1/2,1/2]) spectral gap. 


Proof of Proposition \8.15 Let C\, denote the constant in Propositio n |8.10[ Note 
that C\, was chosen to be sufficiently large in the proof of Proposition |8.10| and can 
be arranged to be taken so that C\, > 9q, where 6*0 is defined by 9q = 2|Vi^i|/g^. 
Also, choose a constant /o such that (q = |Vip|/ 2 g^. Note /o < 9o; below we shall 
see why we make these choices. For (A, (5) in the region (|8.51 ) we have: 


(8.52) 

— det M(£, (5, A, 0) = 7r(£, A)-I-o ((A^-I-£)(A^-I-(5^)) = 7r(£, (5^, A)-b o (£^), 
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where 

7^(£,(5^A)= (g2A2+£Fi,i) (^q^X^ + S^ |Wo,i + hFi,o - . 

We now show that there exists A^’^ S (Co-\/e, ^oV^) such that Tr{e,S‘^, X’^’^) = 0 . 
Note first that eVi_i < 0, £ and the choice of Co implies, upon evaluation of 

7 r(£,( 5 ^, A) at A = Cov^j that: 


^(£, 5 ", Co Vi) = (VCoV + £Vi,i) (g^CoV + <52 < 0 

and by ( 8.52[ ) — det M(£, 5, CoVi, 0) < 0. On the other hand, the choice of 9o 
implies, upon evaluation at A = 0Q^/£ that: 


7r(£, (5^, Oo'/s) — 1 ^ 1 , 1 ) (V^o^ + I WC),i + bbi,o ~ Wi^i | ^ > 0 

and hence, by ( |8.52[ ) — det M(£, 5, OqVs, 0) > 0. Now det M(£, S'^, A, 0) is, for all 0 < 
£ < £ 1 , a continuous function of A. Hence, there exists A®’'^ G (CoVi; ^oV^) such that 
det M{e,S, X^’^, fj, = 0) = 0. Hence, E‘^’^{X‘^’^) = G ifc||= 27 i -/3 ~ spec(-H('^’^l). 
This completes the proof of Proposition |8.15| □ 


Appendix A. Evaluation of £Hi.i for two examples 

Recall the bases {vi, V 2 } of Ah = Zvi 0 Zv 2 and {ki, k 2 } of A^ = Zki 0 Zk 2 , in¬ 
troduced in SectionMore generally, introduce a lattice spacing parameter, a > 
0, and define the scaled lattices: and with bases: v“ = a(^, 5 )^, vf = 

, -i)^ and k? = f , k^ = Uh ^^ere g = ^ and = 

2TrSij. Now introduce the base points: = (0,0) and = a(;^,0) and 

the honeycomb structure with general lattice spacing parameter, a > 0: = 

(A(“1 0 aJj“^) U (B(“1 + a|)”^). To be consistent with previous notation, we write: 
Ad) = A = (0,0), Bd) = B = (^, 0), A^ = A, (A^'V = (A/.)*, Hd) = H. 

Let go(x) denote a smooth, real-valued, radially symmetric (go(x) = go(|x|)) and 
rapidly decaying function on Below we shall use the 2D Poisson Summation 
formula: 

(271)2 


E 


/(x 0 nv) = 


|D/, 




mk-: 


We next present two examples: sums of translates of go over the scaled triangular 
lattice, a|j“\ and honeycomb structure, Hd). In both cases, is expressible in 

terms of % Therefore, if goiO changes sign, then the sign of can be 

changed by varying the lattice constant, a. 


A.l. Example 1: Evaluation of £Vi.i for V equal to a sum of translates 
over the scaled triangular lattice, Define H(x;a) = X)v 6 a'“' (x 0 v). 

The potential H(x;a) is a honeycomb potential; it is periodic, inversion 

symmetric and TZ— invariant with respect to the origin of coordinates Xq = 0. 


Claim 1. HVi = 


(a) _ ( 271 -)^ 


IWI 


ffo 


' ki+k 2 '' 
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Proof of Claim 1. The Poisson summation gives 

( 271)2 


P(x;a) = 


|f2/, 


^ go(mk(“))e* 


(mk(“>) 


Recall also that mk(“) = mik^“^ +TO 2 k 2 “^ = +m 2 k 2 ). Claim 1 now follows 

from: 

y(x;a) = ^ ^ ^ ^ 


lez^ 


V 


A.2. Evaluation of eVi i for V equal to a sum of translates over the 
scaled honeycomb strnctnre H As earlier, take = (0,0)^ and = 
a(^,0)^. The point at the center of hexagon, immediately northeast of is 
Define 

V (x; a) = ^ 5 o (x - A^^^ + ro^“^ + v^ + ^ go (x - + v^ . 


veA 


(a) 


veA 


(a) 


P(x; a) is a honeycomb potential; it is A/j— periodic, inversion symmetric and 7Z— 
invariant with respect to the origin coordinates, Xq = 0, located at the center of a 
hexagon. 

Claim 2. = xso(lii)- 

Proof of Claim 2. Poisson summation yields 


y(x; a) = ^ [gb (■ - -L (mkC))+fo (• - bO + ro^C'j (mkO)] ehmk<“))-x^ 


|!^hl 


riez2 


Now, go (• - A(“) + To(“)) (?) = exp (z? • to(“)) go(?) and go (• - B(“) + (?) = 

exp (i? • (—B(“^ + ro^“^)) (?o(?), and therefore 


P(x;a) = ^^ ^ Jg 2 mk(“).ro<“> ^g 2 mk(“).(-B(“)+ro<“))j ^^(mkW) 




\n 


rieZ2 


Noting that mk^“) • to^“^ = mk • tq and mk^“^ • (—B^“^ + to^“^) = mk • (—B + tq) 
(independent of the lattice constant, a), we have 

^ 27r 

mk^“^ • = (miki + TO 2 k 2 ) • tq = —(2toi - m 2 ), and 

o 

-♦ Ott 

mk^“^ • (-B*^“^ + = (miki + TO 2 k 2 ) • (-B + tq) = - —(2mi - m 2 ). 

o 

Recall again that mk(“^ = mikj“^ + = ^(miki + m 2 k 2 ). Hence, 

(27r)2 ;V +"^2k2 


P(x;a) = 


|H/, 

(27r 


St 


9o 


V 




|H/, 


rieZ2 


2 cos f ^(2mi — m 2 ) 


90 


miki + m2k2 


imk 


(a).. 


imk 


= E e- 

meZ2 

Therefore, 2cos (f) x go {^) 


- ^ 90 


{ 4tt 1 '] 
\V3a) ■ 
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